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A STUDY OF THE SPECTRA OF B STARS 
By OTTO STRUVE 


ABSTRACT 


Thirty-five spectrograms of thirteen stars, ranging in spectral type from Og to B8, 
were measured for the purpose of obtaining a reasonably complete /ist of absorption lines 
in the hottest stars. The list contains 379 lines (between \ 3820 and A 4924), of which 
about 150 are new. Approximately 130 lines have been newly identified. The following 
elements are present: H, Het, Hen, Cu, Cm, Nu, Nim, Ou, On, (Nai), Mgu, 
Alu, Alm, Siu, Sim, St31v, Pu, Pm, Su, Sm, Ceu, Tin, Cru, Fen, Nau; 
and the following are suspected: A 1, Ca m1, S Iv. 

The contours of the Het lines were investigated. Members of the series 2P—mD 
and 2p3— md; were found to be highly susceptible to broadening by electric fields of ions 
and electrons in the reversing layers (Stark effect). Accurate wave-lengths of twenty 
Het lines are given. The wave-lengths of two forbidden lines were found to agree well- 
with laboratory determinations. 

The intensities of several pairs of lines of He t were estimated in 116 stars of types O to 
Bo. The ratio of intensity singlet /triplet is greatest for stars of type B2 and decreases 
rapidly toward class O and slowly toward Bo. This peculiar behavior of the singlets is 
not due to temperature-ionization, nor to effects of pressure or of electrical fields. It is 
possible that the effect is caused by a change of the atomic absorption coefficient depend- 
ing upon outside conditions, or by a departure from thermo-dynamic equilibrium. 

The members of the Fowler series of He 11 were found to be broadened by Stark effect, 
in agreement with theory. 

Several lines of V u and of O 1 which are diffuse in the laboratory (except under 
very low pressures) show no evidence of broadening in stellar spectra. 

The members of doublets and of quartets in O 11 show approximately the same rela- 
tive intensities in the stars as they do in the laboratory. 

The lines of silicon were used for a test of Milne’s generalized ionization equations. It 
is shown that lines belonging to two stages of ionization of silicon are simultaneously 
strengthened in giants. The earlier form of the ionization theory fails to explain this 
phenomenon. It is probably a consequence of Milne’s so-called ‘“‘null-effect.” 

An attempt was made to estimate an upper value for possible deviations in the abun- 
— of silicon in individual stars. The ratio does not exceed 5:1, and probably remains 

elow 2:1. 


It is the purpose of this investigation to provide a foundation for 
further work on the spectra of early-type stars. A preliminary ex- 
amination of spectrograms of several B-type stars showed that there 
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were many absorption lines which did not appear in any of the pub- 
lished lists of stellar wave-lengths.' I therefore decided to measure a 
number of spectrograms in order to have a reasonably complete list 
of absorption lines present in classes O and B. 

The observational material consisted of twelve plates taken with 
the Coudé spectrograph of the Mount Wilson Observatory, and of 
twenty-three plates taken with the Bruce spectrograph of the Yerkes 
Observatory. For nearly all of the spectrograms I used Eastman 
Process plates, which give excellent contrast and bring out many 
lines which are invisible on fast photographic emulsions such as 
Eastman-4o. 

Table I gives a list of thirteen stars selected for measurement. It 
was my intention to cover, as nearly as possible, all subdivisions of 
the Harvard classification between Oe5 (or Og) and B8. The selec- 
tion was based almost exclusively upon the sharpness of the lines. 
Not many years ago such a selection would have been regarded as 
not sufficiently representative of the whole spectral class, because 
there are many stars which have wide and diffuse lines. It is now 
regarded as almost certain that the dish-shaped contours of these 
wide and diffuse lines are produced by rapid axial rotation—a factor 
quite extraneous to the structure of the reversing layer. It is there- 
fore safe to say that our list is representative for normal spectra of 
class B. It includes a number of stars of high luminosity, such as 67 
Ophiuchi and 6 Orionis, and others of comparatively low luminosity, 
such as y Pegasi and ¢ Herculis. 

The plates were measured in the usual way and reduced by means 
of Hartmann formulae. I found it expedient to measure the wave- 
lengths of the faintest lines in one direction only. If all lines had 
been of the same intensity, this procedure would have involved no 
systematic errors of any sort, because the personal difference in 
setting upon stellar absorption lines and upon emission lines of the 
comparison spectrum is eliminated in applying the correction for 
radial velocity. But the lines are not all of the same intensity. There- 

«F, E. Baxandall, Comparison of the Spectra of Rigelian, Crucian and Alnitamian 
Stars (Solar Physics Committee), 1914; H. H. Plaskett, Publications of the Dominion 
Astrophysical Observatory (Victoria), 1, 351, 1922; F. Henroteau and J. P. Henderson, 


Publications of the Dominion Observatory (Ottawa), 5, 5, 1920; B. H. Dawson, Publica- 
tions of the Astronomical Observatory of the University of Michigan, 2, 159, 1916; etc. 


A STUDY OF THE SPECTRA OF B STARS 227 


fore, the possibility remains that faint lines are affected by an error 
of setting: the procedure is strictly correct only for lines of the same 


TABLE I 
OBSERVATIONAL MATERIAL 
Linear 
Sp. oe a Measured | No. of 
Star Plate Date UT. Range Meas- 
44500 (AA) ures 
Oo IR 1930 Aug. 23 3hg5™] Y.O. 40”! 3958-4713 I 
IR 9506 Sept. 4 840 | ¥.0. «40! 30 3958-4713 I 
9 «Cas... Bo 1B 1727* | 1908 Sept. 8} 15 10 | ¥.O. 40”! 30 3999-4572 I 
1B 2482* | 1910 Aug. 12 | 19 46 | Y.O. 40} 30 4572-4982 I 
37 Ori... Bo IR 9661 1931 Feb. 25 3990-4172 I 
237 Sco... Bo 1R 9709 Apr. 11 g2t | Y¥.O 30 4012-4713 I 
28 CMa.. Br R 1771 1930 Mar. 3 rar | ¥.0; 4338-A7590 I 
1R 9307 Mar. 6 040 | Y.O. 40”! 30 3901-4717 I 
8 BCep... Br IR 9548 Oct. 19 | 3032-4713 I 
88 y Peg... Ba B 246* | 1901 Nov. 27 | 12 11 | Y.O. 40”! 10 4172-4885 I 
tR 9249 | 1929 July 14 | 700] 40”! 30 3013-4164 | 1 
IR 9253 July 17 | 30 3761-4274 2 
R 1588 Sept. 7 408 | Y.O. 40”! 10 4338-4713 1f 
R 1599 Sept. 21 438 | 10 4338-4722 2 
C 342 1930 Oct. 1 909 | Mt.W.100”} 12 3840-4603 I 
C 345 Oct. 2 745 | Mt.W.100”) 12 3840-4716 2 
826 Cet... B2 1B 4323* | rors Dec. 14 | 1316 | Y.O. 40”! 30 3958-4555 I 
IR 9513 1930 Sept. 9 Sar ac") 30 3058-4713 I 
24 y Ori.... B2 C 63 1927 Jan. 23 8 00 | Mt. W. 100” 3 4338-4495 I 
C ay 1929 Oct. 23 | 11 03 | Mt. W. 100” 6 3957-4603 2 
Cas... B3 IR 9269 July 29 | 811 | Y.O. 30 3901-4384 1** 
85. Her.... B3 IR 9242 July 3 36 3913-4172 2t 
© aie 1930 June 18 755 | Mt.W.100”} 12 3930-4405 I 
C 314 June 19 5 27 | Mt. W.100”| 12 3930-4495 I 
325 July 5 740 | Mt.W.100”) 12 3903-4713 4§ 
67 Oph.... Bsp IR 9256 1929 July 19 1 30 3013-4481 I 
IR 9262 July 27 3013-4481 2t 
1930 June 19 | 9 43 | Mt.W.100”| 12 3930-4405 I 
C 326 July 6 653 | Mt. W.100”| 12 3903-4727 4 
19 B Ori.... B8p C 122 1928 Dec. 12 7 33 | Mt. W. 100” 6 4005-4737 I 
C 129 Dec. 30 7 50 | Mt. W. 100” 6 4326-4057- I 
C 362 1929 Mar. 26 | 317 | Mt. W. 100” 6 4320-4057 I 
R 1768 1930 Feb. 20 006 | Y¥.O. I0 4338-47590 I 
1R gs51s Sept. 9 | 10 40 | Y.O. 40”! 30 3913-4717 I 
IR 9516 Sept. 9 | 1054 | Y.O. 40”! 30 3958-4716 I 


* The time is G.M.T. counted from noon. 

+ Measured by F. E. Roach. 

t One measure by P. C. Keenan. 

§ One measure by F. E. Roach. 

** Only a few of the stronger lines measured on this plate. 


intensity as those from which the correction for radial velocity was 
derived. But in practice this does not appreciably affect the results. 
The probable error of setting on a faint line so much exceeds any 
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systematic error that the latter may be disregarded. The wave- 
lengths of the stronger lines were determined in the usual way from 
measurements made in two directions, and should be accurate within 
a few hundredths of an Angstrom unit. 

The wave-lengths of the lines and their estimated intensities ap- 
pear in Table II. Intensity 1 denotes a very faint line which can be 
seen only with difficulty. A faint, but easily visible, line is denoted 
by 3. The strongest lines, other than those of hydrogen, have an 
intensity of approximately 1o, and the hydrogen lines are estimated 
at about 20. Considerable attention has been devoted to the faintest 
lines: those of intensity 1 were carefully scrutinized and rejected 
from the table unless they appeared on several plates or unless they 
were definitely identified with elements already known to be present. 

The work was confined to the ordinary photographic region of the 
spectrum, which was further limited on the violet side by the trans- 
mission curve of the spectrograph. The list of wave-lengths is reason- 
ably complete between \ 3950 and \ 4700. The lack of sensitivity 
of the Eastman Process emulsion to radiations of wave-lengths 
longer than \ 4700 made it necessary to employ. a coarser emulsion 
for the region between d 4700 and ) 4900; and consequently only 
the stronger lines could be measured. In the case of y Pegasi the 
measurements extend on the violet side to \ 3850. 

Of the thirteen stars in Table I, ten are separately listed in 
Table II. The remaining three stars—¢ Cassiopeiae, ¢' Orionis, and 
y Orionis—were used in determining the mean wave-lengths but 
are not given in the section entitled “Intensities.”’ g' Orionis was 
measured in detail over a narrow range of wave-lengths, and in 
¢ Cassiopeiae and y Orionis no attempt was made to measure the 
faintest lines. 

The wave-lengths in Table II are the arithmetical means of all 
independent determinations. Unit weight was assigned to each meas- 
urement, irrespective of the intensity of the line. This procedure was 
adopted to avoid throwing too much weight upon any one meas- 
urement. Some of the best plates available, principally those from 
Mount Wilson, were measured several times. Since each measure- 
ment was given weight 1, this somewhat increases the effect of the 
better plates upon the adopted wave-lengths. 
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TABLE II 


WaveE-LENGTHS OF LINES IN B-TypE STARS 


INTENSITY 
AVE- se a a] 
3820. . 15 Het 9.61(4); Het 9.75(1) 
20 Hn 5.39 
3850.94 I Ou 1.04(3); Ou o.81(2); On 
1. 58(2) 
3856.01 2 Ou 6.16(5); Nm 6.07(3); Sim 
6.02(8) 
3857.10... I Om 7.18(4) 
3862.53... 2 Sim 2.59(6) 
3864.45... I Ou 4.45(5); Om 4.68(1) 
3867.55... 9 Het 7.46(2); Het 7.62(1) 
3869 .88.. I [Al 1 0.06(3)] 
3871 .88.. 2 Het 1.80(1) 
I 4. 10(2) 
6.41(2); Cm 6.67(1); Cu 
6.05(1) 
3889.02 9,055 Het 8.65(10); Ci 
9.18(1 
3809 . 53 [Al 11 0.68(10)] 
3907 . 49. 2 r |. O 1 7.45(4) 
3911.99. 7 3 |...]...].-- |---| Om 1.95(10); Om 2.09(2) 
3918.95. 6 4 414 |...| Cu 8.98(6); Om 9.28(6); Nu 
9 .00(6 
3920.65.. 5 514 |...| Cm 0.68(8) 
3923. 30.. I a 1 Sm 3.47(4) 
3924.42.. 4 4.44(4) 
3920.59.. 4 313 |---| Her6.53(2) 
3928.59.. 2 Su 8.59(5) 
3933 -60.. 8 | 6 |10 | Cam 3.66(10) 
3930.05.. | 2] Hers.o1(1) 
3938. 20.. 1 |...| [Alm 8.62(})] 
3045.01.. O 1 5.05(5) 
3954. 30. . 4/3 | 2]. O 11 4.37(7) 
3955 .80.. 2 N 11 5.85(6) 
3961.52 3 ...| 1.59(8); S mt 1. 56(3) 
3964.67...| 2 815 |6|4|6|7|4 | He14.73(4) 
3970.16...|20 |...|...130 |40 |23 |...|15 |22 |20 | He 0.08; S 11 0.72(3) 
$14] 2 & | Omg. 


Ce 
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TABLE II—Continued 


INTENSITY 
OBSERVED 

3979.67... | Str 9.85(2) 
3980. 24.. Al 1 0.56(2); C 0. 35(2) 
3083 .73.. I z (Al 3. 
3985 .80.. 2 I S mt 5.95(4) 
3003. 33-- I 2 | 2 | S11 3.49(6) 
3994.098...| 2 | Nm §.0o(10); [Alm 5.86(s)] 
4009.28...} 3 | Herg.27(r); 9.90(2)] 
4024.05. I Het 3.99(1); Om 4.04(1) 
4026.21.../10 | 8 {15 |12 |15 |15 | 8 | Het 6.19(5); Het 6.36(1); Alu 

6.5(5); 6.00(3) 

4028.73 1 | 2 | S11 8.80(5) 
4030.51 I 1 | | Almo0.87(0); Alm 1.14(4) 
4035.12...]...|..., 2] N115.09(4); O11 5.00(0) 
1] 1] g.30(4); Al 119.40 (0) 
21 |] On 1. 20f0) 
4043.46 
4048 .15 I 2 I O 1 8. 22(1) 
4057.09.. (1) 1] 1 | Put 7.39(4); Nm 7.00(1) 
4000.74...| I O 11 0.58(3); O11 0.98(2) 
4062 .90.. O It 2.90(1) 
4064. 89.. I I & 4.40(3) 
4067 .04.. I Nim 7.04(3) 


On g.90(6); Om 9.64(4); Cur 
0.30(3); Cm 8.97(3); Su 
9.65(1); S 9.89(1) 
4072.11...) Om 2.16(8); Pi 2.13(5) 

1] 1 Om §.87(10) 
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TABLE Il—Continued 


INTENSITY 
ales 
aie 
4101.76...|20 | 8 |30 |20 |20 |30 |50 |50 |30 |20 
4120.84...) 6} 81 e138] 6) 3 
I 


4143.80... 


(3) 


IDENTIFICATION 


..| 8.86(4); O 11 g.00(0) 


P 0.04(7) 


.| , 


| N 1 2.28(2) 
O 11 3.90(2) 


..| Owms.1a(3) 


| Om 7.16(2); N 7.35(0) 


Si Iv 8.86(10) 


| Pm 1.53(s) 


..| O11 g.28(4); Si tv 8.86(10) 

| 

| 


| Om 2.94(5); [N 1 2.65(0)] 

[O 11 5 .63(0)] 

Ou 7.25(4); Nim 7.31(10); 
O 7.32(1) 

S 1 9.42(3); S 1 9. 27(2) 

He 11 0.08 

Hé 1.74 

N mt 

O11 5.00(7); O11 4.74(5) 

Ou 7.07(1) 

O 110. 80(3); 0. 20(1) 

2.04(4) 

Si tv 6.10(8) 

O 11 9. 22(8) 

Het 0.81(3); Het 0.98(1); Ou 
0.30(3); On 1.48(4); OM 
0. 55(2) 


| Siu 8.05(8); P 1 7.49(6) 


| Ou 9.34(2) 


Al 0. 21(2) 


..| (0110. 74(0)] 


| Het 3.77(2); Ow 3.77(2); On 


| 
| Sr 2.31(4) 
| 
| 3.52(1); forbidden Het 3.55; 


5.90(0) 


0.88(10); S11.0.90(2) 
Er | | 4 
2|...| Sm5.16(4) 
rj...]...]...| Om 6.09(3); Nu 5.76(3); On 
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TABLE II—Continued 


INTENSITY 
OBSERVED 
LA 

a a 2) is] ve) ~ 

Al 9.92(1) 


ARG 1 |...| Om 6.54(3); [N  6.8(0)] 
2] 2/213 ]2 | Sm 2.73(6); Sa 2.37(3) 


106.00... 
4169.04...| 2 |X | 2 


312]3 | 2] Om 3.31(7); Su 3.14(6) 
[N 11 5.0(0 


4|2|2| Het 8.97(1); Om 9.23(4); Alu 
8.51(1); Al 8.42(3) 
4.07(4) 

2 ; ...J...| Om §.45(8) 


I 
2 

2 


274...) 2 1(3))-2 | ]...| O11 g.79(10); S 11 0.07(2) 

4t09.92...| 2 |(2)] 2 Hem 9.87; m1 0.02(6) 


~ 


~ 
a 
= 
= 
mn 
= 
mn 
= 


6. 35(0) 

2] 7.83(3); Abi 6.82(6); 
Alt 7.49(5); Al 7.98(4) 

4237.01...|...]...) 2] 2 2 | N 11 6.98(6) 


4241. (1)| N 1 1.80(8) 


an 
w 
~ 


| 
| 
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TABLE IIl—Continued 


OBSERVED 
WAVE- 
LENGTH 
LA. 


INTENSITY 
Qa 
gic 9 
el al Bilal Ss] 
2| 2 
I 2 
aha 
2 
I 
I 
Sear ae te oe 
2 I 
.|20 |10 |20 |20 {15 |40 |40 |50 |24 |20 
2/3/13 


IDENTIFICATION 


Ci 7.27(10); Ci 7.02(8); Su 
7.86(5); Sm 7.2744); Su 
6.90(4) 

S g.80(2) 

[O 1 3.17(0)] 

O11 5.52(4) 

O 11 6.71(1) 

7.40(1); O11 7.90(1) 

S 1 8.62(3) 

O11 1.40(0) 

Om 2.96(1); Su 2.65(3); Ou 
3.75(0); Al 2.97(3) 

Sm 5.00(5); 5.70(3) 

Ou 8.83(1) 

Om 1.25(1); Om 2.23(0) 

S 4.44(6) 

O 11 4.82(3) 

Fe 11 6.56(6) 

[Tz 11 0.05(12)] 

Feu 3.18(4); S m 3.78(1) 

O 1 3.82(5) 

O 1 7.31(1); [Al 1 7.20(3)] 

Ti 7.89(12) 

O 1 2.10(0) 

O 3.43(1) 

[Ti 1 4.98(5)] 

O11 5.55(0) 

O 1 7.16(8) 

Su 8.66(3) 

O 11 9.65(8); O 11 9.93(1) 

O11 §.77(3) 

O 1 7.48(3); 7.89(0) 

Om 1.89(2); Ou 1.47(0); Ou 
1.21(0) 

S ut 2.69(6); Om 2.76(1) 

O 11 6.86(6) 

Ti 11 7.92(10) 

Hy 0.47; Ou 0.36(2); Sur 
0. 25(6) 

O 11 2.00(4); O 1 3.30(0) 

Om 5.57(7) 


| 
| 
4253.66. Ou 3.98(8); S mt 3.51(6) 
4207.10. 
4269. 43. | 
4271.78. | 
4273.35. 
4275.68.. 
4276.95. | 
4277.55. 
i 4278.60. 
4281.22. 
4283 .04.. 
| 4285.14. 
4286. 24.. 
4288 .93.. 
4291.48. 
4294. 27.. 
4294.87. 
4200. 54.. 
4208.94. 
4299.88. 
4303. 28.. 
4303.90. 
4307 .37- 
4308.42. 
4310. 54.. 
4352.22. 
4313.33-- 
4314.44.. 
4315.27. 
4317.21... 
4318.45... 
4319.66.. 
4325.64.. 
4327.75.. 
4331.79.. 
4332.68. . 
4330.91... 
4337.85... 
434°.47.. 
4342 .90.. 
4345. 50.. | 
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TABLE II—Continued 
INTENSITY 
OBSERVED 
Blo 
~|/Ola 
4347-47 I 2 I Om 7.43(5); Alt 7.78(4); Alu 
7 .80(3) 
2] 1] Omg.44(8) 
4351.27 2 | 3 | 2 Om 1.28(6) 
4358.06...|...]...]...] 2 O 7.25(0); O11 8. 40(0) 
4361.51...) 2] 1% S ut 1.57(6) 
4372.30 Om 1.65(2); Cu 2.49(1); Cu 
1. 59(0) 
4379.22 4 2 N ut 9.09(10); O 9.55(3) 
I .|...| [Forbidden He 1 3. 26] 
4387.93-..|9 | 6| 8 lio 17 |10 15 | 9 | 4 | Het 7.93(3) 
4390.55 2 2/1 | 2] Mgtro.58(10) 
113 ]2] 2] 2] 2 O115.95(7) 
4405.99 I O 11 6.02(1) 
4409 . 32 Cr 1.00(a); Cm 1.5a(2); Cu 
0.06(1) 
4414.90 Om 4.80(10); [Pu 4. 28(6)] 
4416.74 ...| 2 | Fem 6.81(4) 
4416.97 1 | 2 |(2)) 6.97(8); [Pu 7. 30(5)] 
ape Nu 2.71(6); Sm 2.45(2); [Alu 
| | 2.82(4)] 


| 
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TABLE IIl—Continued 


| 82 Cet Bo 


| 85 Her B3 


INTENSITY 

E- Nn 

4437-55---| 2} 6 
4471.56...|11 |10 |10 |10 | 8 |17 


BASS 


4489.51... 


4401 


2 2 I 
I 


4534-99... 


| 67 Oph Bsp 


12 


| 


| 19 B Ori B8p 


IDENTIFICATION 


N u 3.48(2); Mg 1 3.90(8) 

Het 7.55(1) 

N ur 1.99(3) 

O 3.05(5) 

Nu 7.04(10); Om 7.08(1); 
[Al 1 7.8(3)] 

O 8. 21(6) 

O 2.38(6) 

Sm 3.66(3); Sm 4.45(3) 

Ou 5.40(4); [V 1 5.54(0)] 

O 11 6.32(2) 

Ou 7.88(4); Om 9.32(3); Sur 
7.76(1) 

Forbidden He I 9.97 

Het 1.48(6); Het 1.69(1) 

P 11 5.26(8) 

Ou 7.88(2); Nu 7.74(2); Sm 
8.44(1) 

Al ut 9.97(4); Al 1m 9.89(3) 

Mgu 1.33; Mgm1.13 

S m1 3.48(4) 

Ou 8.09(2); O11 9.48(1) 

Fe it 9.21(4) 

O 11 9.48(1); [Al 1 9.87(3)] 

Ou 1.25(3); Al 0.g0(2) 

Feu 1.41(4) 

Ti 11 1. 27(25) 

O 6.50(2); N 7.58(3) 

N 7.58(3); Om 6.50(2) 

Fe 8. 29(8) 

N ut 0.92(6) 

Al 1m 2.54(4) 

N ut 4.80(7) 

Fe tt 5.33(6) 

N 8. 18(3) 

Fe 11 0. 24(6) 

Fe  2.64(6) 

N 1 3.60(4) 

S 1 4.89(6) 


Al int 9.18(6); Al mt 8.91(1) 
N 11 0.37(5); Pu 0. 78(7) 
Mg 1 4.26(4) 

N ut 4.57(3); S Iv 5.60(4) 


| 
gar 
I 
12 |r0 | 
I 
| 8 | 
la | 
Leu 
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TABLE IIl—Continued 


INTENSITY 
OBSERVED 
LA. 

Blo 

ASAT 03 1 | I He 1.63 
4549.44. ..|(1)|. 1 | 3 | 3 | Femg.48(4); Su 9.54(4) 
4552.04 jas 5 |2| Simt 2.65(9); Nim 2.50(4); 

Sm 2.44(5) 


Si 1 7.87(7) 


4588.15 ..| | Alt 8.19(5); Al 8.08(3) 
4590.97..-.|1 14/3 | 3] 2 O110.98() 

4590.14. O 11 6.19(8) 

4609 70 > Om 9.42(4); Nu 9.60(0); Om 


0. 14(3) 


4613.78. 

4619.15. 

4629.44... | Fem 9.33(4); P 8.70(5) 
4630.55...|...| 3 3 5 | 3 |---| 4 


4634.50 2 2 I I .| N mt 4.16(8) 
4638.71 1}/4/2/4/3]4] 2 O 11 8.86(6) 
4640.57 2 Alt 0.36(4); 0.38(3); 


N ut 0.64(10) 
1 ]...]...] Om 1.83(9); Nim 0.64(10) 
N ut 1.90(3) 
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TABLE II—Continued 


INTENSITY 
OBSERVED 
3 g IDENTIFICATION 
LA. 
2 
2|1 | Omg.15(10) 
Ou 0.85(6); Alm 0.54(2); Al 
0. 65(r) 
4683. 76. I Si 3.77(2) 
4685.88...} 8 5 | 1 Het 5.81 
1.76(0) 
4702.092...)...|.-.| 1} 2 I O 11 3.18(3); N 11 2.57(0) 
4705.36...| 1 |(2)] 2 | 3 3 Om §.36(8) 
|10| 7 |X |X | 4] Her 3.14(3); Her 3.37(1) 
2 | 6.66(s); Sm 6. 26(4) 
NOTES 
3928.59 If 67 Ophiuchi is included, the wave-length is \ 3928.46. 
3933.14. This line is blended with Ca 11 \ 3933 on most plates. 
3933.60 Wave-length based upon 8 Orionis alone, since in other stars it is of inter- 
stellar origin. 
3966.91 Measured twice on the same plate. Uncertain. 
3968.46 Wave-length based upon 8 Orionis alone, since in other stars it is of inter- 


stellar origin. 
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3980. 
3979- 
3989. 
3991. 
4015. 
4030. 
4045. 
4953. 
4064. 
4067. 
4069. 
4084. 
4097. 


4099 . 
4124. 
4129. 
4129. 
4135. 
4140. 


4140. 
4150. 
4153. 


4159. 
4164. 


4186. 


4188. 
4199. 
4202. 
4207. 
4220. 
4227. 
4232. 
4243. 
4263. 
4275. 
4276. 
4277. 


27 
15 
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Probably blended with S 1 \ 3979.85(2). 

Ingram (Physical Review, 32, 172, 1928) gives intensity 4. 

Uncertain. 

Uncertain. 

Attributed to Ni mu by Baxandall. 

Several close A/ 1 lines. , 

There is a forbidden line of He 1 at \ 4045.17, but the identity is improbable. 
There is a forbidden line of He 1 at \ 4053.56, but the identity is improbable. 
Uncertain. 

Attributed to Ni 1 by Baxandall. 

10 Lacertae gives \ 4070.10 and may be a blend with some other lines. 
Perhaps a blend. 

N ut probably predominates in 10 Lacertae, while O 11 predominates in 
6 Canis Majoris. 

Uncertain. Very close to //6. 

Uncertain. 

Probably partly blended with Si 11 \ 4128.05. 

Uncertain. 

Uncertain. 

Rufus, Sawyer, and Paton (Publications of the Astronomical Observatory of the 
University of Michigan, 3, 261, 1923) mention a line of Si at \ 4141.04. 
Ishida gives a forbidden line of Het at X 4141.32. 

Uncertain. 

Uncertain. Perhaps A/ ur \ 4149.90 (3), etc. 

The wave-length is very sensitive to spectral type, because O u predomi- 
nates in the early stars, while S 1 is more prominent in the later ones. The 
following wave-lengths are obtained separately for each subdivision: 


4153.23 


Uncertain. Several other faint lines of A/ u in this region. 

The wave-length does not agree with that of the S 1 line. The intensities 
suggest that the stellar line is mostly due to some other element. 

The wave-length in y Pegasi is 6.44, agreeing roughly with K 1 6.24(8), the 
strongest line of ionized potassium. The wave-length in 1o Lacertae is 6.96. 
Uncertain. 

The stellar line is probably a blend. 

Uncertain. 

Uncertain. 

Uncertain. Baxandall gives a line at \ 4220 in ¢ Orionis. 

Several other faint lines of A/ 1 in this region. 

Uncertain. 

Uncertain. 

Very doubtful. There is a line of K 1 at 3.40(7). 

There is a forbidden line of He 1 at 5.81, but the identity is improbable. 
Uncertain. 

The line of 6 Ceti is probably blended with S 1 8.51(4). 
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84 
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95 
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4278. 
4283. 


42806. 
4294. 
42096. 
4299. 


4303. 
4308. 
4314. 
4325. 
4327. 
4332. 
4347 - 
4350. 
4351. 
4352. 
4358. 
437°. 
4395. 


4400. 
4404. 
4405. 
4409. 
4447. 
44606. 
4401. 
4527. 
4529. 
4549. 
4553- 
4572. 
4589. 
4597. 
4602. 
4604. 
4011. 
4619. 
4630. 


4034. 
4638. 
4649. 


4058 
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The line of y Pegasi is probably blended with O 1 7.40(1). 

In 8 Canis Majoris and 8 Cephei O 1 predominates, the stellar wave-length 
being 3.29. In 67 Ophiuchi and c Herculis S m predominates, and the stellar 
wave-length is 2.82; y Pegasi gives an average wave-length of 3.01. 

This line seems to be blended in y Pegasi with \ 4185.14. 

The wave-length in y Pegasi is 4.47, the line being blended with O 11 4.27. 
The line in y Pegasi and 7 Scorpii is probably not due to Fe 11. 

The individual wave-lengths are very discordant, and there is some doubt as 
to the reality of this line. 

Possibly blended with O 11 3.82(5). 

Uncertain. Baxandall gives in 8 Orionis 8.1. 

Very uncertain. 

The line in B Orionis, 5.84, is probably not due to O 1. 

Very doubtful. 

Uncertain. 

Several other faint lines of A/ 1 in this region. 

Somewhat uncertain. 

The line in 67 Ophiuchi, 1.50, is blended with Fe 1 1.77(10). 

Blended with 3.84. 

Uncertain. 

Uncertain. 

In B Orionis, 5.43, and in 67 Ophiuchi, 5.76, this line is blended with 77 u 
5.04(25). 

The identification is doubtful. 

Uncertain. 

Uncertain. Perhaps identical with O 1 6.02(1). 

Uncertain. 

This line is blended with O 11 8.21(6) on Yerkes single-prism spectrograms. 
In 6 Canis Majoris this line is blended with O 11 5.40(4). 

This line is partly blended with O 1 9.48(1). 

Uncertain. 

On some of the one-prism Yerkes plates this line is blended with NV 11 0.37(5). 
The wave-length in 10 Lacertae is 9.11; this line is not due to Fe 11. 
Uncertain. 

The identification is uncertain. 

Very doubtful. 

Uncertain. 

This line is blended with N 11 1.49(8). 

Uncertain. 

Uncertain. 

Rufus, Sawyer, and Paton (op. cit., 3, 261, 1923) give a line of Si at \ 4619.60. 
In 7 Scorpii and 8 Cephei this line is blended with Si iv 1.38(3); and in 
B Orionis, with Fe 11 9.33(4). 

Partly blended with unknown line at 3.75. 

Rufus, Sawyer, and Paton (0. cit., 3, 261, 1923) give a line of Si at 4638.36. 
This line is blended on some of the plates with C 11 0.24(2). 

Uncertain. 
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4666.12 Very doubtful. Rufus, Sawyer, and Paton (oP. cit., 3, 261, 1923) give a line 
of Si at \ 4665.76. 

4678.28 Very doubtful. 

4709.54 Uncertain. Rufus, Sawyer, and Paton give a line of Si at \ 4709. 20. 

4726.59 Somewhat uncertain. 


No attempt has been made in Table II to provide accurate stellar 
wave-lengths which would replace those now being used in determi- 
nations of radial velocities. The purpose of this table is to serve as a 
basis for identifications. A comparison of my stellar wave-lengths 
with values obtained in the laboratory shows that the agreement is in 
general very satisfactory, the difference Star — Laboratory rarely ex- 
ceeding a few hundredths of an Angstrom unit. While the few lines 
showing greater discordances may be due in part to errors of meas- 
urement—either in the stellar spectra or in the laboratory plates— 
some of the discrepancies are doubtless caused by erroneous identi- 
fications or by the blending of adjacent lines on spectrograms of 
small dispersion. 

The probable identifications are given in the last column of 
Table II. In the case of blends, the components are named in their 
probable order of importance. Laboratory lines in square brackets 
indicate improbable identifications. Additional information not con- 
tained in the body of the table may be found in the notes following it. 
If a line measured in a given star was subsequently regarded as 
suspicious, the estimate of its intensity was placed in parenthesis 
and its wave-length was not used in forming the mean. A line which 
was definitely seen, but not measured, is indicated by a cross (X) in 
the list of intensities. 

The total number of lines in Table II is 379. Of these, approxi- 
mately 150 are new. Nearly 130 lines have been newly identified 
with known elements. The number of unidentified lines is 73. 

Table III contains a list of all elements found in stars of class B 
and the number of lines identified for each element. Group I gives 
the number of lines which are either free of blends or which appear as 
stronger components in blends with other lines; group II gives the 
number of lines which occur as fainter components in blends with 
other lines of the same ion; and group III gives the number of lines 
which are of uncertain identification or which are blended with 
stronger lines of other elements. 
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TABLE III 
NUMBER OF LINES 
ELEMENT Source FoR LABORATORY-LINES 
I Il Ill 
1 Giiiicneateseoe A. Fowler, Report on Series in Line-S pectra, p. go. 
21 6 3 | Lbid., p. 93. 
gl 5 Renner: 8 7 7 | A. Fowler and E. W. H. Selwyn, Proc. Roy. Soc., A, 120, 
312, 1928. 
Cur 2 I o | I. S. Bowen and R. A. Millikan, Phys. Rev., 26, 310, 
1925. 
RES Sa8is.0 aN 3 15 | A. Fowler and L. J. Freeman, Proc. Roy. Soc., A, 114, 
662, 1927 
4 | L. J. Freeman, ibid., A, 121, 318, 1928. 
3 eee IOI 25 21 | A. Fowler, ibid., A, 110, 476, 1926. 
O 1 | Ibid., A, 117, 317, 1928. 
iS) 2 I A. Fowler, Series in Line-S pectra, p. 120. 
Alu ’, 6 | 21 | R.A. Sawyer and F. Paschen, Annalen der Physik, 84, 
I, 1927. 
ALE 2b. 8 4 3 | F. Paschen, zbid., 71, 142, 537, 1923. 
1 | A. Fowler, Phil. Trans. Roy. Soc. of London, A, 225, 
I, 1925. 
Soils Tbid. 
g | I. S. Bowen, Phys. Rev., 29, 510, 1927. 
Pim C8 A ete R. A. Millikan and I. S. Bowen, ibid., 25, 600, 1925; 
M. O. Saltmarsh, Phil. Mag., 47, 874, 1924. 
20 4} 19 | M. J. Gilles, Recherches sur la structure du Soufre 
(Paris, 1930), S. B. Ingram, Phys. Rev., 32, 173, 1928. 
Pb re WO 5 | M. J. Gilles, of. cit., S. B. Ingram, op. cit. 
1 | M. J. Gilles, op. cit. 
1 | T. L. de Bruin, Proc. Roy. Acad. (Amsterdam), 29, 713, 
1926. 
A. Fowler: Series in Line-S pectra, p. 127. 
2 | J. A. Anderson, Ap.J., 59, 76, 1924; I. S. Bowen, Phys. 
Rev., 31, 497, 1928. 
Tiu 2 | H.N. Russell, Ap.J., 61, 223, 1925. 
T. Dunham, Jr., The Spectrum of a Persei (Princeton, 
New Jersey, 1920). 
Feu 7 H. N. Russell, A p.J., 61, 223, 1925. 
CINE A. G. Shenstone, Phys. Rev., 30, 255, 1927. 


HYDROGEN 


The hydrogen lines have been discussed at some length in a recent 
paper by Elvey and Struve." The most interesting results are the 
relation between line-contour and Stark effect (illustrated in Figures 
8 and g of that paper), and the pronounced tendency for stars of 
high absolute magnitude to have Balmer lines of comparatively small 
total intensity. The contours of Hé in 67 Ophiuchi (giant) and 85 c 


* Astrophysical Journal, 72, 277, 1930. 
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Herculis (dwarf) are reproduced in Figure 1. 67 Ophiuchi is of a later 
spectral type than c Herculis, and we should expect the hydrogen 
line to be stronger in the later type, but the opposite is observed. 
The same effect may be seen in the halftone illustrations (Plates 
XI-XV). The Balmer lines are especially sharp and narrow in 
8 Orionis, while in y Pegasi and in c Herculis these lines have ex- 
tended wings. 
NEUTRAL HELIUM 

A p pearance of lines—Neutral helium is represented in the B-type 

stars by the following lines: 


He LINEs IN B STARS 


Triplets Singlets 
DIFFUSE SERIES DIFFUSE SERIES 
PRINCIPAL SERIES PRINCIPAL SERIES . 
SHARP SERIES SHARP SERIES 


I have shown in several articles' that the members of the diffuse 
series of both singlet and triplet systems are broadened by the action 
of ionic fields in the reversing layers. The amount of broadening 
depends upon the density of the charged particles—ions and free 
electrons—and is therefore a measure of the star’s luminosity. The 
evidence can now be extended to higher members of the series. In 
the reproduction of the spectrum of y Pegasi, Plate XI, the lines 


' [bid., 69, 173, 1929; 70, 85, 1929; 70, 237, 1920. 
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NA 3820, 3872, 3927, 4009, and 4026 appear conspicuously broad 
and diffuse, while AX 3867, 3878, 3935, and 3965 are much narrower. 
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Fic. 1.—Contours of Hé6 in 67 Ophiuchi (dots) and 85 « Herculis (circles). The 
ordinates are percentages of absorption and the abscissae are Angstrom units. 
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Fic. 2.—Microphotometric tracings of spectra of 67 Ophiuchi (top) and of 85 ¢ 
Herculis (bottom) made with the Coudé spectrograph of the Mount Wilson Observa- 


tory. 


‘Microphotometric tracings of three plates are shown in Figures 2 
and 3. The curve for y Pegasi shows conspicuous widening in \ 3927 
and \ 3872, while \ 3867 remains narrow and deep. In ¢ Herculis, 
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d 3927 is very broad, and resembles the same line in y Pegasi. On 
the other hand, 67 Ophiuchi, a giant, shows a narrow line at \ 3927, 
indicating little, if any, Stark effect. 

The contours of the helium lines are shown in Figure 4 with two 
other lines for comparison. The relation of line-contour to spectral 
series is again well marked. 

The broadening of the diffuse members can be ascertained also in 
d 4388 and \ 4472, as I have already shown in earlier papers; and in 
agreement with J. Pauwen' the forbidden line of He, \ 4470, is seen 


Hn H¢ He 
| | | 


He 3820 

Te 3868— 
Te 3872— 
2027 —. 
Te 3936— 
He 3065 
He 4009 — 
He 4026— 


Fic. 3.—Microphotometric tracing of spectrum of y Pegasi taken by Dunham 
with the Coudé spectrograph of the Mount Wilson Observatory. 


to be broader than lines of other elements in the same region of the 
spectrum. 

Wave-lengths of the He lines.—Table IV contains the wave-lengths 
of the helium lines as measured on the best available plates. The 
values agree, in general, with those in Table I, which are based upon 
more extensive material. 

In a previous paper? I made the suggestion that the wave-lengths 
of the helium lines should be affected by the character of their Stark 


- patterns. Certain lines are shifted by electrical fields toward the red, 


while others are not appreciably displaced, and a few are shifted 
toward the violet. It should be emphasized that the effect, if actual- 
ly present, would be very small. As was shown in the case of the 
hydrogen lines,3 the upper layers of the star’s atmosphere super- 


1 Ibid., 70, 263, 1929. 2 [bid., 70, 92, 1929. 
3 Elvey and Struve, Astrophysical Journal, 72, 277, 1930. 


| 
| 
| 


A STUDY OF THE SPECTRA OF B STARS 245 


impose upon the broadened wings a sharp and narrow core that is 


not affected by electrical fields. This core would, therefore, not be 


shifted at all; and it is only in the wings that a lack of symmetry 
might be detected. 

The wave-lengths in Table IV indicate several cases of abnormally 
large differences Star—Laboratory. In the stars the line d 3964 is 


rest T re T T T T T T 

e 


Fic. 4.—Contours of lines in y Pegasi: 
He 4009 He 3927 He 3872 
He 4026 He 3820 He 3965 
He 3872 Ci 3919 3921 
Abscissae are Angstrom units; ordinates are percentages of absorption. 


systematically shifted toward the violet, while \ 4169 is shifted 
toward the red. This is in excellent agreement with my former re- 
sults based upon single-prism Yerkes spectrograms.’ However, I am 
not now certain that these shifts are actually due to Stark effect. 
The line \ 3964 falls upon the wing of the hydrogen line e, as may 
be seen on the microphotometer tracings in Figures 2 and 3. It is to 


* [bid., 70, 94, 1929 (Table I). 
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be expected that a narrow line falling upon the wing of a broader 
one, would be shifted in wave-length away from the center of the 
broader line. The reason for this is the increased hydrogen absorp- 
tion on the red side of the center of the helium line as compared to 


TABLE IV 


WaAvE-LENGTHS OF He I 


| 
| 
= = = = S 
3867 .46\ (2) | 
.351| .347| .237| .369] .321]..... 
4026. 189\(5) | 
4o26.358/(1) .242| .241| .230] .216| .240]..... 
120.812 | 
4387.928 (3) .952 935) .9140.997, .978.0.941 0.963 
4437-549 (1) .640} .580] .583] .581]..... .498| .601] .650)..... 
4469.97 .834| .g80] .898] .964]..... | 
.559| ©.440,0.450 .595| -544| .634 
13.143\ 
4921.929 (4) | | 0.899 0.028 
| | 


that of the violet side. I do not know whether this effect would ex- 
plain the observational data quantitatively, but it may, at least in 
part, account for the shift. The line \ 4168.97 is a blend, and is 
therefore not suitable for a test of the shift resulting from Stark 
effect. The oxygen line \ 4169.23 has a laboratory-intensity of 4 and 
would tend to shift the blended line in the direction observed. While 
it is probable that some of the stellar wave-lengths of He 1 actually 
differ from those established in the laboratory, we shall leave the 
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question open as to whether or not appreciable shifts are caused by 
Stark effect. 

Forbidden lines of He 1.—-The strongest forbidden line of He 1, 
d 4470 (2p3—4f3), has been measured in many stars, and a reliable 
stellar wave-length is therefore available for comparison with the 


laboratory results. 
STELLAR WAVE-LENGTH 


The whole material.................... 4469 95 


The best laboratory determinations, as quoted by Y. Ishida and 
T. Tamura,’ are: 
LABORATORY WAVE-LENGTH 


4469.97 
4469.99 


The agreement is excellent, and there can be no doubt as to the 
correctness of my identification. 

In addition, one other forbidden He 1 line has now been identified 
with certainty. Its stellar wave-length is 4025.32 (2p'— 5f5), agreeing 
well with the following laboratory values: 


sfs 583 
Ishida and Tamura........ 4025.58 Ishida and Tamura........ 4025.48 


The list of forbidden lines observed in the laboratory by Ishida 
and Tamura gives the approximate electric field at which each line 
first becomes visible. The lowest value, 100 volts, is assigned to the 
line \ 4470. Its strength in the stars is thus accounted for. There is 
another laboratory line, \ 4007.79 (2P—7P), which appears in a 
field of 140 volts. I have measured on several plates a line of mean 
wave-length \ 4007.53. The agreement is satisfactory, but the stellar 
line is so faint that its identity cannot be regarded as definitely 
established. It is much fainter than the other forbidden lines, 
d 4470 and \ 4025, the latter of which is first seen in a field of <400 
volts. At my request Dr. Ishida has checked his observations of 

* Scientific Papers of the Institute of Physical and Chemical Research, 13, 118, 1930. 
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d 4007.79, and he finds that the limiting field strength, which was 
intended to give only a rough indication of the conditions, may be 
somewhat larger than 140 volts. However, he is convinced that it 
should not be higher than 400 volts/cm. Professor T. Takamine has 
also informed me that the line at \ 4007 is easily produced in the 
laboratory. More extensive data will be needed in order to decide 
whether the discrepancy between laboratory and stellar intensities 
is real. 

Of the other forbidden lines listed by Ishida and Tamura, none 
would be expected to occur in stellar spectra. Nearly all of the 
stronger forbidden lines are so close to the permitted lines of He 1 
that they are blended. A few doubtful coincidences in wave-lengths 
are recorded in Table II. 

Intensities of He lines.—The relative intensities of the helium 
lines are not the same in all stars. This was first noticed some years 
ago’ when I compared the estimated intensities of the triplet line 
d 4472 with those of the singlet line \ 4388. I have now extended 
the material to include other lines, omitting all stars showing rapid 
rotation. Unfortunately, the number of stars with perfectly sharp 
lines is small, and it was therefore necessary to use a number of stars 
that had moderate rotational broadening. 

The following results are intended to serve as a basis for a more 
precise spectro-photometric study of the helium intensities. I have 
therefore made estimates of the relative intensities of the following 
three pairs of lines: 


4388 (2P—sD) _ 4009 (2P—7D) (2p'— 5s!) 
4472 (2p'—4d3) 4026 (2p'— sd’)? 4144 (2P—6D) 


All stars have been reclassified according to the Harvard (and 
partly Mount Wilson) criteria. It will be seen that the agreement 
with the spectral types given in the Henry Draper Catalogue is (with 
a few exceptions) very good. I have attempted, furthermore, to 
make a rough classification within each spectral subdivision, ac- 
cording to absolute magnitude. This division into three groups— 


™ Nature, 122, 994, 1928. 
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giants, intermediate stars, and dwarfs—is admittedly uncertain, and 
depends primarily upon the following criteria: 

1. Presence or absence of wings in the Balmer lines (cf. y Pegasi— 
dwarf; 8 Orionis—giant). 

2. Presence or absence of forbidden line of He 1 at \ 4470. 

3. Appearance of helium lines belonging to diffuse series. 

4. Intensity of interstellar lines of Ca 11. 

Unfortunately, all four criteria are affected by rotational broaden- 
ing, and this causes considerable uncertainty in the estimates. The 
latter are shown in Table V. The rotational broadening has been 
estimated on an arbitrary scale, zero indicating perfectly sharp lines 
and to standing for appreciable broadening. It should be remem- 
bered that stars with rapid rotation were not used in the work. The 
relative intensities of the lines are designated as follows: 

= equal 
< slightly fainter 
<< fainter 


<<< much fainter 
<<<< very much fainter 


An inspection of the table gives the following results: 

1. In spectral subdivisions Br, B2, and B3, all stars resemble 
rather closely the spectrum of y Pegasi (Plate XI). The singlet lines’ 
are here nearly as strong as the triplet lines. It will be recalled that 
in the laboratory the singlets are usually much weaker than the 
triplets. There seems to be no difference, in these spectral classes, 
between giants and dwaris. 

2. As we proceed to the earlier spectral classes, Bo and O, the 
singlet lines fade out much more rapidly than the triplet lines. There 
are indications that this is more pronounced in the giants than in 
the dwarfs. 

3. Among the later spectral classes, B7, B8, and Bo, there is a 
similar fading-out of the singlets with respect to the triplets, but the 
effect is not as pronounced as in classes O and Bo. Supergiants of 
type B8, such as 6 Orionis, show less fading-out of the singlets than 
do stars of intermediate or low luminosity. 

The halftone reproductions of stellar spectra in Plates XI-XV illus- 
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TABLE V 


EstTiIMATED RELATIVE INTENSITIES OF HELIUM LINES 


Star | Mg. | Harvard 4388 | 4009 
| | Sp. | | 4472 | 4026 4144 
Spectrum O 
Giants: 
e208. ....: | + 9°52’| 3.7 | Oes | 2] <<< | <<<< > 
a2 2 +61 48 | 5.2 Oes5 | 4 << K<<< = 
Intermediate: 
omon.....| 6 35 + 9-50 Oes |4| << Rik 
Dwarfs: 
7 CMa... .. —24 46 | 4.4 | Oes5 | 6 << 
TO 22 26 +38 32 | 4.9 | Oe5 | < << < 
Spectrum Bo 
| 
Giants: | 
oem. +66 10 | 4.4 Bo 4 
Intermediate: | | | 
27 +62 23 | 4.2 Bo | 3 < 
+34° 980 Aur.| 5 10 | +34 12] 5.8] Bop | 1 < <<< | = 
— 116] 1.8 Bo 5 < 
AB 5 34 — 2 39 | 3.8| Bo << | 
53 x Ori. 5 43 — 9 42] 2.2] Bo 
Dwarfs: | 
27 — 7 23 | 4.6 B3 ° < << | = 
16 30 | —28 2.9] Bo o| = << | = 
59° 2334 Cep..| 21 9 | +50 34] 5.6] B2 | 3] < << | = 
Spectrum Br 
Giants: | | 
3 48 +31 35 | 2.9 | Br 3| = < 
47 p-1L00......| 10 28 + 9 49 | 3.8| Bop | 3 = |< < 
Intermediate: 
5 29 + 9 25 | 4.5 Bo o| = = 
Bond 708..... 5 30 — § 290 | 6.3 Bo 4 = << = 
26 EMa.:...| 6 18 —17 54| 2.0] Br I = < < 
| (aoe 6 49 —20 6| 4.7 Br 3 = < < 
—28 50 | 1.6 Br I = < 
16 15 —25 21 | 3.1 Br 5 = < 
Dwarfs: 
—13 17] 4.3 | Br I os = 
—1°935Ori...| 5 28 — 140| 5.3] Ba o| = < < 
—6°1234..... 5 30 — 6 4] 4.7 Br o| = << > 
4#CMa..... 6 28 —23 21 | 4.4 Br 2 > < = 
Ce... +70 7 | 3.3] Br o| = < < 


Star 


Giants: 
62 x7 Ori..... 


Intermediate: 


Per... 


..... 


55 Cyg....... 
114 Tau...... 
—6° 
—1° 1004 Ori. . 


—19°1407CMa 
20:4 LyF........ 
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TABLE V—Continued 
Spectrum B2 
| 
| | +20° 8’| 4.7| Bap | 1 = < < 
21 35 +61 38 | 4.9 Bap | 2 = | < < 
| 336 | +33 30|5.0] Bz |4/ = |< = 
8 40 —32 50 | 3-7 B2 c = |< < 
17 16 | —24 54 | 3.4| 1| = l< < 
22 37 | +39 42 | 5.2| Ba | = | < < 
| o 8 +14 38 | 2.9 B2 | ° = | < < 
| O31 | +53 21 | 3.7| B3 o| = | < < 
| 2 34 6| 4.0 B2 ° = < 
4 31 — 3 33 | 4.1 B2 o| = | < < 
| 4 46 + 5 26] 3.8 B3 2 = | < < 
| § £7 — 0 29 | 4.6 B3 ° = < < 
| + 616] 1.7 B2 3| = l< < 
| 18 4 +20 48 | 4.3 B3 2| = < < 
Spectrum B3 
20 46 +45 46 | 4.9 B2 | 1 = < < 
238 | +2717 | 4.6 B3 7 1 < < 
449 | +217|3.9] B3 = < < 
Sad | | 4.8} Be = < < 
5 26 Bz | 4 = 
5-30 | — 6. I = < < 
5 36 +16 29 | 4.9 B3 8 < << < 
5 36 —11r|5.0| B3 5 = < < 
6 14 —19 56 | 5.3 B3 3 = < | < 
Ig 10 +38 58 | 4.5 B3 = < < 
Ig 28 +34 14 | 4.8 B3 2 = < < 
19 34 + 5 10] 5.2 B3 9| = < < 
20 27 | +48 37| B3 | 7] = < < 
20 39 | +490 5.4 B3 2/ = < < 
2148 | +48 51/4.3| B3 | 3] = < < 
22-17 +1I 42 | 4.9 B3p | 5 = < < 
22 26 | +42 37|4.5| B3 | 4] = < < 
I 47 +63 II | 3.4 B3 ° = < a 
a Rte +20 47 | 5.2 B3 2 < << = 
| 3 40 + 5 44| 5.4 B3 6 < << < 
35% | +s 5-3); BS 3i = < < 
413 | +5041 | 5.5| B3 | 4] = < < 
457 — 719 4.8 B8 2 < < < 
| 457 | +58 50/5.3| Bap | 6| = |< < 
| 59 +41 6 | 3.3 B3 7 < << < 
5 2 | +4 Bs | 3) = < < 


42 0 Oph. . 
Dwarfs: 
88 y Peg... 
$2.6 Cet. ... 
yp ..:.. 
230 
24.5 On... 
In 
| 
440 Aql..... 
45 w' Cyg.... 
Stet 
an Pee. 
Dwarfs: 
46° 5. 
+50° 973 Per 
65.0 Er... .. 
Cam... ..... 
107 Aur..... 
Taw... 
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TABLE V—Continued 
Star | a | 5 Mg. Harvard Rot. 4388 | 4208 452i 
Sp. 4472 4026 4144 
Spectrum B3—Continued 
Dwarfs:—Cont. 
Siep........| | —24° | Bs I = < 
67:2) 6 2 +14 47 | 4.4 B2 = < < 
—6° 1446 Mon| 6 7 — 6 32] 5.1 B3 o| = < < 
re CMa...... 6 16 —30 1 |.3.1 B3 2 = << < 
10 Mon... 6 23 — 4 42| 5.0 B3 4 = < =< 
17 58 +20 50] 5.1 B3 o| = < < 
58 | 36 sx B3 5 = < 
+50° 2708 Dra! 18 58 +50 23] 5.2 B3 3 = < < 
Ig 20 | +29 26] 4.9] B3 = < < 
Cyg...... 21 39 +50 44 | 4.9 B3 4| = << 
FO UAE. 22:52 +41 4] 5.5 B3 ° = < < 
Spectrum B4, Bs, Bo 
Giants 
25 x Aur 5 26 +32 71] 4.9 Br ° < < = 
20+ CMa..... 6 52 —16 55 | 4.4 | Bs I = < 
2407 CMa....| 6 59 —23 41 | 3.1 Bsp | 1 < << < 
17°50 | 2:56 | | sp | x = < < 
Intermediate 
— 6 57 | 3.7 Bs << <<< = 
5 34 | +25 50| 5.0] B3 4| = < 
—4° 1362 Mon| 6 2 —A tr | Be < << < 
BOMEX. IO 25 —o 7|5.0| Bs << << 
Dwarfs: 
29 m And..... © 32 +33 10] 4.4] B3 o| = eg > 
AT... ... 244 | 3 Bs < 
BY. 3 49 —24 54 | 4.8 Bs Al <<< 
—16° 796 Eri 4 —16 39 | 5.4 | B3 << = 
eS 4 14 +46 16 | 4.9 B3 ° < KK = 
20 p Aur +41 42] 5.1 B3 5 = 
eg48U....:.. 5 42 +13 52] 5.2 Bs 4| < << < 
—8° 4010..... 15 20 — @ sr | 5.2) << << > 
— 227 Her 16 17 +46 33 | 3.9 Bs O} 1< << > 
— 85. Her +46 4] 3.8 B3 o| = < 
3 Vul 19 19 | +26 41] 4.9] Bs5 a} = << 
20 28 | +10 58| 4.0] a <<< = 
30 Peg 22°15 +517] 5.4 Bs 3 << < 
2 Lac 22.17 +46 2] 4.7 Bs AP) << << < 
Spectrum B7, B8 
Giants 
zo. | 8 19 | BSp < < = 
18 8 —ar ‘< < = 
+46° 3111 Cyg] 20 52 | +47 2] 5.9] B8p | 1 = << > 


| | } | | | | | 
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TABLE V—Continued 


Star a | Mg. | Harvard 4388 
| 4472 4026 4144 
Spectrum B7, B8—Continued 
| 
+53° 102 Cas.| ob3r™ | +53°37’| 5.1 Bs < << 
46 w Cas...... 1 48 +68 5.0 B8 3 << 
526 17 | 4.4 Bs 2 < <4 
114-6: 26 +28 31 | 1.8 B8& 3 < 
+16° 794 Tau.| 5 26 | +16 59] 5.5 Bo a) “<< << 
+65 50 | 3.2 Bs << < 
+10 55 5.1 B8 << 
Spectrum Bo 
Giants: 
2 Hev Cam:...| 3) 21 +59 36 | 4.4 | Bop ° < << 
22 20 | +48 58 | 4.6 | B&p | o | < | <<< > 


trate these points. Consider first the spectrum of y Pegasi (Plate XI), 
a typical dwarf of class Bz. Here \ 4388 is approximately equal in 
intensity to \ 4472, and X 4oog is slightly fainter than \ 4026; simi- 
larly, \ 4121 is slightly fainter than \ 4144. 

Consider next the two stars « Herculis, a dwarf of class B4, and 67 
Ophiuchi, a giant of class Bs (Plate XII). The relative intensities of 
d 4009 and \ 4026 are the same in both stars, and are about equal 
to those of y Pegasi. We conclude that there is no effect of absolute 
magnitude in the intermediate spectral types. 

On the other hand, in 10 Lacertae, a dwarf of type Og (Plate XIID), 
d 4009 is extremely faint, while \ 4026 is fairly strong. However, 
since 10 Lacertae is characterized by much greater ionization than 
y Pegasi, all its neutral helium lines are comparatively weak, which 
suggests that the relative intensity of singlets and triplets may be 
merely a function of the total intensity of the helium lines: 


Ising] 
=f (Itriplet) 
[triplet 


That this is not the whole truth may be ascertained from the spec- 
trum of 6 Orionis (Plate XV), a supergiant of class B8. Here \ 4026 
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is weaker than it is in 10 Lacertae, but the intensity of \ 4009 is 
much greater. We conclude that while the ratio of the singlets to 
the triplets decreases on both sides of the maximum of the helium 
lines, this decrease is much more rapid on the side of high tempera- 
tures than on that of low temperatures. 

It is difficult to explain this peculiar behavior of the helium lines. 
It is not due to Stark effect, since it does not depend upon absolute 
magnitude in the intermediate spectral types B1 to Bs. For the 
same reason, pressure is excluded as a possible cause. It is probable 
that the peculiar intensities of the helium lines in stellar spectra are 
of a totally different nature than the changes in the relative inten- 
sities of singlets and triplets frequently observed in the laboratory" 
and invariably found to be related to the gas pressure. 

The intensity of a stellar absorption line depends upon three fac- 
tors: (1) the number of atoms in the upper and lower states which 
define the transition corresponding to the line, (2) the atomic ab- 
sorption coefficient, and (3) the process by which the radiation is 
transmitted through the stellar atmosphere. 

The number of atoms in any energy-level is given by the theory 
of ionization, provided we may consider the reversing layer as being 
in a state of thermodynamic equilibrium. This assumption has 
proved satisfactory in many applications of the theory, and we shall 
adopt it here, subject to further revision. We may then regard the 
number of atoms in the various states as defined uniquely by the 
temperature, the pressure, the potentials of ionization and of excita- 
tion, and by the statistical weights of the states. To be specific, 
consider the two lines \ 4026 (2p’—5d5) and \ 4009 (2P—7D). The 
relative number of atoms in the two states 2p’ and 2P is: 


— 11,600 
T 


'L. S. Ornstein, H. C. Burger, and W. Kapuscinsky, Zeitschrift fiir Physik, 51, 34, 
1928; W. Hanle, ibid., 56, 94, 1929; W. Elenbaas, ibid., 59, 289, 1930; W. C. Michels, 
Physical Review, 36, 1362, 1930. 
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where w is the statistical weight of the state, E the excitation poten- 
tial in volts, and 7 the temperature in degrees Centigrade. For the 
two lines under consideration, 


E,= 20.81 volts , 


E,= 21.12 volts . 


For a temperature of T = 20,000° C, corresponding roughly to spec- 
tral type B2, we have 


The numerical factor differs very little from unity. Even for a low 
temperature of only 10,000° C, it is equal to 1.4. We are therefore 
justified in concluding that the ratio V,/N, is practically independ- 
ent of the temperature and that we cannot account for the observed 
changes in the relative intensities of the two lines by attributing 
them to corresponding changes in V,/N,. 

The process by which radiation is transmitted through the atmos- 
phere of a star does not provide an explanation of the observations. 
Given the atomic absorption coefficient o, the intensity in a line is 
simply 

I=1,f(o) , 


where the form of f is independent of outside conditions, such as 
temperature or pressure. 

There remains the possibility that the atomic absorption coeffi- 
cients corresponding to the two lines might vary with the conditions.' 


* A. Unséld has confirmed, in the case of hydrogen (Zeitschrift fiir Astrophysik, 2, 
199, 1931), the theory of V. Weisskopf and E. Wigner (Zeitschrift fiir Physik, 63, 54, 
1930; 65, 18, 1930), according to which the absorption coefficient produced by radiation 
damping involves not only the oscillatory strength of the line but also the lives of the 
atom in its upper and lower states. In a recent paper Unsild (Zeitschrift fiir Astro- 
physik, 3, 94, 1931) has called attention to the fact that in helium the singlets have a 
very strong transition probability from state 2'p to state 1's, while there is no corre- 
sponding probability for the triplets. Consequently, the functional relation of the ab- 
sorption coefficient to wave-length is not the same for triplets and singlets, the latter 
being relatively the broader. This should, of course, result in a strengthening of the 
singlet lines, provided that all contours are widened predominantly by radiation damp- 
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If this should not prove to be the case, we should be tempted to 
doubt the assumption of thermodynamic equilibrium. 


IONIZED HELIUM 


Ionized helium is represented in Table II by the following lines: 
d 4685.81 (Fowler series); and AA 4541.63, 4199.87, 4100.08 (Picker- 
ing series). The lines \ 4859.36, near HB, and d 4338.71, near Ly, 
were observed in several stars by H. H. Plaskett' but could not be 
seen on the Yerkes spectrograms because the latter were not suffi- 
ciently exposed to bring out faint lines in the immediate vicinity of 
the Balmer lines. 

The two members of the Pickering series, \ 4542 and X 4200, are 
appreciably broadened in many stars. The spectrum of to Lacertae, 
Plate XIII, illustrates this. On the other hand, the Fowler line \ 4686 
is not appreciably broadened. We shall show that the broadening of 
the Fowler lines is caused by ionic Stark effect. 

The He 11 ion is similar in structure to the atom of hydrogen, and 
its two series are expressed by the following relations: 


Fowler series: = 
ow i v=4R 


Pickering series: v= 
v=4R 4 [n'P 


The Fowler series originates from transitions between an upper 


energy-level of principal quantum number n’=4, 5, 6,....anda 
lower level, »=3; for the Pickering series the principal quantum 
numbers of the two levels are n’=5, 6, 7,....and n=4. In the 


presence of an electric field each of these states is split up into a 
number of separate levels characterized by the parabolic quantum 


ing. Unséld suggests (private communication) that the observed differences in the 
ratios, singlets/triplets, may be due to the fact that radiation damping is not equally 
effective in all stars. Lines widened by thermal Doppler effect or by Stark effect would 
not be subject to the influence of the transition probabilities from the lower level. This 
question will be investigated further. 

t Publications of the Dominion Astrophysical Observatory (Victoria), 1, 325, 
1922. 
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numbers k, and k, and by the equatorial quantum number m. Ac- 
cording to quantum mechanics,’ 


n=k,+k.+m+1 
Hence, 
o<ki<n—-1; o<k,<n—-1. 


The displacements of the components in the Stark-effect pattern 
are given by the following formula,’ in the usual notations: 


[n’(ki— —n(k.—k;)] . 


For hydrogen the formula remains the same, except that the coef- 
ficient outside the brackets is twice as large. 
The factor 
3hF 
1677em, 


is independent of the quantum numbers. Consequently, in compar- 
ing the relative displacements of corresponding components of two 
lines of He 01, it is sufficient to consider the expression 


X=n'(ki—ki)—n(k.—k,) . 


The maximum value of this expression for \ 4686 is obtained by 
substituting 


n' = 4; n= 3; ki=3; ki=o; k,=0; k,=2. 


We then obtain 
Xmax(X 4686) = 18 . 


For \A 4542 and 4200 the quantum numbers for maximum separa- 


tion are: 
n k’ ka ki 


4542....8 4 7 ° ° 3 
4200...10 4 9 ° ° 3 


! E. U. Condon and P. M. Morse, Quantum Mechanics, p. 82, 1929. 
2 R. Minkowski, Handbuch der Physik, 21, 407, 1920. 


| 


258 OTTO STRUVE 


Consequently 
X max(A 4542) = 68 


4200) =102. 


If we assume, for simplicity, that the distribution of intensities 
among the various components of the Stark pattern is the same for 
all three lines, we obtain the following approximate ratios of the 
widths of the two Pickering lines to that of the Fowler line: 


Width (A 4542) _ 3. Width (A 4200) _ 


Width ( 4686) Width 4686) 


This presupposes that the widths of all three lines are caused by 


Stark-broadening. 
It is of interest to compare these values with those for the hydro- 


gen lines: 


Ha HB Hy Hd He 

3 4 5 6 7 
2 2 2 2 2 
2 2 3 4 5 6 

I I I I I 
3 Xmax 4 7 II 16 22 


The values of $X max Should be directly comparable to those of He 11. 
We see that the Stark-broadening of both Fowler lines exceeds that 
of He, and would be equaled only by the higher members of the 
Balmer series, where blending and confluence of the wings make it 
impossible to study the contours in detail. 

It may seem surprising that the total width of the hydrogen lines 
greatly exceeds that of the Pickering lines. This is due to the greater 
abundance of hydrogen in the reversing layer of the star. The maxi- 
mum separation in the Stark effect, computed above, affects the 
atomic absorption coefficient and causes the latter to be much 
flatter and broader in the case of He 11 than in the case of H. But 
the contour of a line depends upon the optical depth, and conse- 
quently involves the number of atoms per cubic centimeter. There- 
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fore, the wide wings of the less abundant element He 11 may easily 
appear fainter than the wings of the more abundant H. 

The appearance of the lines of He 11 in 10 Lacertae agrees quali- 
tatively with the computed values. This star has probably a com- 
paratively dense atmosphere, since the forbidden line of He 1, \ 4470, 
is present and the wings of the hydrogen lines are broad and strong. 
It would be of interest to investigate the appearance of the Fowler 
lines in stars of low density. Unfortunately, there is appreciable ro- 
tational broadening in all the brighter stars of class O which, from 
the appearance of the hydrogen lines and from the strength of the 
interstellar line of Ca 11, are suspected to be supergiants. In the 
star g Camelopardalis, however, where the rotational velocity is not 
large, the Fowler lines seem to be appreciably narrower than in 10 
Lacertae. 

An effect of broadening similar to that observed in 10 Lacertae 
was noticed by F. Paschen' in the laboratory for the higher mem- 
bers of the Fowler series, in the extreme ultra-violet part of the 
spectrum (AA 2733, 2511, 2385, 2306, 2253, 2215, 2186, 2165). The 
broadening of the Pickering lines \ 4542 and \ 4200 has not, to 
my knowledge, been observed in the laboratory. 

The relative intensities of the Fowler line \ 4686 to those of the. 
Pickering lines seem not to be the same in all stars. While in 10 La- 
certae \ 4686 is very strong, it is faint in g Camelopardalis. 


SINGLY IONIZED NITROGEN 


There are 49 lines of N 1 in Table II. Of these, 31 are free from 
blends or constitute the stronger components of blends. It is well 
known that many lines of N 1 are diffuse in the laboratory,? and 
W. E. Pretty? has shown that in some of the lines this broadening 
is accompanied by a shift in wave-length. The whole phenomenon 
may, according to Pretty, be due to Stark effect. F. E. Baxandall* 
believed that those lines which are hazy in the laboratory have also 


' Sitzungsberichte der Preussischen Akademie der Wissenschaften, p. 135, 1926. 

2 A. Fowler, Proceedings of the Royal Society of London, A, 114, 664, 1927. 

3 Proceedings of the Physical Society of London, 41, 442, 1929. 

4 Researches on the Chemical Origin of Various Lines in Solar and Stellar Spectra 
(Solar Physics Committee), p. 52, 1910. 
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a tendency to be hazy in stellar spectra, and I have suspected this 
also." 

Table VI shows the lines of intensity 3 or greater given by A. Fow- 
ler as being diffuse in the laboratory. A careful examination of these 
lines on my best plates fails to show any evidence of broadening 
(see y Pegasi, Plate XI). Baxandall mentions also’ the line \ 4530 as: 
being broadened in stellar spectra. I believe that this may now be 
definitely ascribed to blending with Al 1, 4529.09. The presence of 


TABLE VI 


Hazy LIneEs oF N 11 


A Int. Remarks 

3n Blended with He 1 4026. 19(s5) 
4n Blended with O 11 4035.00(o) 
5n Blended with O 11 4041. 20(0) 
3n 

gn 

3n Blended with several lines of 

Alu 

6n Blended with S 4432.45(2) 
3n 

3n Not observed 


the latter line in stars having strong lines of N u (vy Pegasi and prob- 
ably also y Orionis) was not known to Baxandall. 

Of the lines in Table II, there is only one, \ 4227.83, which, ac- 
cording to Pretty, shows an appreciable pressure-shift in wave- 
length. The stellar wave-length agrees exactly with the value found 
in the laboratory under low pressure. It is probable that the electric 
fields in stellar atmospheres are too small to produce measurable 
effects in nitrogen. While no numerical values are available, I be- 


_ lieve that this is what should be expected from the amount of Stark- 


broadening observed in helium and hydrogen. The Stark effect de- 
creases very rapidly with increasing atomic number. 


SINGLY IONIZED OXYGEN 


Singly ionized oxygen is represented in Table II by 147 lines, of 
which 101 are either entirely free from blends or constitute the major 


' Astrophysical Journal, 70, 96, 1929. 2 Tbid. 


| 
| 
| 
if 
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portions of blends, 25 appear as fainter components of blends with 
other lines of Ou, and 21 are present as fainter components of 
blends with lines of other elements or are identified but tentatively. 

The spectrum of O1 has been analyzed by A. Fowler,’ R. H. 
Fowler, and D. R. Hartree;? H. N. Russell;3 I. S. Bowen;? C. Mihul;5 
and F. Croze and C. Mihul.® It consists of doublets and quartets 
and contains numerous intercombination lines. The intensities have 
not been determined accurately in the laboratory, but there are 
available estimates, on an arbitrary scale, by A. Fowler.’ Many lines 
appear broad and diffuse in laboratory spectra (except when the 
pressure is very low). From analogy with .V i, it seems probable 
that Stark effect is responsible for this broadening.* 

The great number of O 1 lines in stellar spectra makes this ele- 
ment particularly suitable for detailed study, and for this purpose 
more accurate determinations of the laboratory intensities will be 
necessary. We shall at this time consider only two problems in con- 
nection with O 11: the relative intensities of quartets, doublets, and 
intercombination lines in the laboratory and in stars; and the ap- 
pearance in stellar spectra of lines which are diffuse in the laboratory. 

We have seen before? that in the case of helium the singlets are, 
as a rule, much stronger in the stars than in the discharge tube. It_ 
is of interest to find whether such a difference exists in the case of 
Ou. Table VII gives a summary of the estimated intensities of lines 
not seriously affected by blends. The difference between doublets 
and quartets is very slight. The number of lines due to intercombi- 
nations of the two systems is not sufficient to permit a reliable con- 
clusion, but it is probable that there are no large differences between 
laboratory sources and stars. 

Table VIII contains all lines of laboratory intensity 3 or greater 
between A 3850 and \ 4700 which, according to A. Fowler, are 
diffuse (marked ‘‘n’’ in his paper). The only lines free from blends 


' Proceedings of the Royal Society of London, A, 110, 476, 1926. 


2 Thid., A, 111, 83, 1926. 5 Comptes rendus, 185, 937, 1927. 
3 Physical Review, 31, 27, 1928. 6 Thid., 185, 702, 1927. 
4 [bid., 29, 231, 1927. 7 Loc. cit. 


8 W. E. Pretty (Joc. cit.) found reasons for attributing this type of broadening in 
N 11 to Stark effect. 
9 Page 240. 
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are AX 4275, 4294, 4303, 4378. Our spectrograms (see especially 
y Pegasi, Plate XI) give no definite evidence that any of these lines 
are broadened in the stars. There is a suspicion of faint wings on 


TABLE VII 


INTENSITIES OF Ou LINES 


vy PEGASI B Canis Majors 
Int. | No. | Int. | No. | Int. No. Int. | No. Int No. Int No. 

TABLE VIII 
Hazy LINEs or O11 
r Int. Remarks 

ESS 3n Blended with O 11 4060. 98(2) 

4n Blended with S7 Iv 4088 .86(10) 

4n Blended with N 111 4097. 31(10); O If 4097. 32(1) 

8n Blended with S m 4253.51(6) 

4n 

3n Blended with S m1 4285 .00(5) 

5n 

4n Blended with Hy 

2.55. 3n 

4401. 25 3n Blended with A/ 111 4490. 90(2) 

4n Blended with N 11 4609 .60(0); O 1 4610. 14(3) 

gn Blended with O 11 4609.42 


either side of \ 4303, but these may be unreal or they may be due 


to other elements. 
SILICON 


The silicon lines are particularly suitable for a test of E. A. Milne’s 
generalized ionization theory,’ because the majority of B-type stars 


' Philosophical Transactions of the Royal Society of London, A, 228, 421, 1929. 
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show two stages of ionization and in some stars even three stages are 
represented. 

Consider the two stars 67 Ophiuchi and « Herculis. The first is a 
giant of Harvard class Bs, while the second is a dwarf of class B3. 
Both stars show the lines of Si 11, AA 4128 and 4131, and of Sz 10, 
A 4553, 45608, 4575. In 67 Ophiuchi all five lines are appreciably 
stronger than in« Herculis. This increase in the strength of lines be- 
longing to two successive stages of ionization is unexpected and leads 
to interesting conclusions. According to the older theory, by Saha, it 
has usually been assumed that an increase in intensity of lines belong- 
ing to a certain stage of ionization can only be produced at the expense 
of some other stage. Now, in the case of 67 Ophiuchi and of « Her- 
culis, we are on the high-temperature side of the maximum of S7 1, 
but on the low-temperature side of the maximum of S7 111. A lower- 
ing in temperature alone would necessarily produce an increase in 
the intensity of Si 1 and a decrease in that of Si m1. However, 67 
Ophiuchi and « Herculis differ not only in temperature but in ab- 
solute magnitude as well. Consider, for simplicity, stages 1, 2, and 3 
of the silicon atom. We have, in the designations of Adams and 


Russell," 


28X107 T? u, —11,600/, 


N Pe Ur T 

T3 U2 —11,600/, 
= 

N, 3.28X10 e 


To these equations is added the basic assumption that V,+N,+N, 
remains constant throughout the sequence of stellar spectra. It is 
clear that an increase in NV, must be compensated by a decrease in 
N, or in N,, or in both. But in spectral classes B3 and Bs the lines 
of neutral silicon have disappeared, and it is difficult to escape the 
conclusion that an increase in NV, must be compensated by a de- 
crease in V,. It would thus appear that the theory of ionization in 
its older form does not account for the simultaneous increase in the 
intensities of Si u and Si 11. 

We have implied in the preceding discussion that the depth of 
the reversing layer is the same in all stars. It remains to be seen 
whether a general increase in optical depth would account for the 


™ Astrophysical Journal, 68, 11, 1928. 
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observations. A general increase in optical depth would enhance all 
lines, including those of hydrogen. As a matter of fact, the hydrogen 
lines decrease markedly in intensity as we pass from ct Herculis to 
67 Ophiuchi. Similarly, the helium lines are appreciably stronger 
in « Herculis. However, it is probable that this is due to Stark- 
broadening. 

The difficulty may be overcome by applying the new theory of 
Milne. For the two stars in question we find surface gravities of ap- 
proximately 103 and 10! cm/sec?. Consequently, Milne’s' computa- 
tions and diagrams for silicon are immediately applicable. Accord- 
ing to his scheme, the increase in intensity of Si 1 would be account- 
ed for by the lower temperature of 67 Ophiuchi. Because of the so- 
called “‘null-effect”’ of Milne, there is no effect of absolute magnitude 
for Si 11. On the other hand, the lines of Si 11 are very sensitive to 
absolute magnitude, and this not only would balance the expected 
decrease in intensity due to temperature but would even result in a 
strengthening of Sz 11 in 67 Ophiuchi as compared to ¢ Herculis. 

The same explanation may be applied to the two stars: y Pegasi, 
a dwarf, and g Cephei, a giant. Both stars are classified as B2. In 
the giant, Siu and S71 are greatly strengthened, and I suspect 
that this is accompanied by a slight strengthening of the lines of 
Sitv. If the temperatures of these two stars were the same, Si II 
should not be strengthened. Since it actually is stronger in g Cephei 
than in y Pegasi, we should conclude, on the basis of Milne’s theory, 
that 9 Cephei has the lower temperature. It is difficult to test this 
in other lines without a detailed knowledge of their relation to 
absolute magnitude. 

If we consider next two stars of type Og to Bo, namely 1o Lacertae 
and g Camelopardalis, the former a dwarf, the latter probably a giant, 
we obtain a similar result for $7 11 and Si 1v. The lines of Sz 111 are of 
about equal intensity in the two stars: Milne’s null-effect on the 
high-temperature side of the maximum of S7 t11 prevents a change of 
intensity depending upon absolute magnitude. For the same stars, 
however, Silv is on the low-temperature side; consequently the 
giant 9 Camelopardalis has much stronger lines, AA 4089 and 41106, 
than 10 Lacertae. 

Loc. cil., 48§2. 
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The foregoing comparisons are purely qualitative and will in due 
course be replaced by actual measurements of line-contours. It is of 
interest, however, that Milne’s generalized theory appears to give a 
much better representation of the observed intensities of the silicon 
lines than does the older form of the ionization equations, which 
dealt exclusively with concentrations without reference to atmos- 
pheric depth.’ 

An interesting side-issue is the determination of a lower limit for 
possible differences in total abundance of silicon in the atmospheres 
of stars of any one spectral subdivision. It is significant that when- 
ever stars of similar spectral type show different intensities of the 
silicon lines, Milne’s theory provides a satisfactory explanation. If 
differences in abundance were large, we should expect that stars, 
assigned to the same spectral class on the basis of their oxygen or 
helium lines, might have totally different lines of silicon. In classify- 
ing the spectra for Table V, I have attempted to use as many criteria 
as possible. The question arises whether these spectra are homo- 
geneous except for effects of absolute magnitude. Let us consider, 
for example, the dwarfs of spectral class B2, listed in Table V. The 
lines Si 1 4128 and \ 4131 are not sensitive to small differences 
in absolute magnitude and are therefore suitable for this test. An 
intercomparison of these lines in all eight stars shows very small 
differences, which may perhaps be explained by small errors in the 
classification. However, to obtain a lower limit, let us suppose that 
all the changes noted are caused by differences in the abundance of 
silicon. In order to obtain numerical values, we note that the dif- 
ferences in the intensities in no case exceed the normal difference 
between the two multiplet lines \ 4553 and \ 4574 in any one star, 
and that in most cases they are well below the difference between 
44553 and 4568. The lines \X\ 4553, 4567, and 4574 form a 
triplet of Si 1m (23s—3%p), and have the following relative inten- 
sities in emission: 5:3:1. These values are proportional to the num- 
bers of atoms contributing to the three component lines.’ It is rea- 
sonable to conclude that individual stars do not differ in abundance 

"See also Publications of the American Astronomical Society (Forty-sixth meeting), 


7, No. 2, 1931. 
2 Struve and Elvey, Astrophysical Journal, 72, 267, 1930. 
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of silicon by more than a ratio of 5:1 and that the ratio is usually 
well below 2:1. 
OTHER ELEMENTS 

Of the other elements identified in B-type stars S 11 leads with 
43 lines. S ut has 15 lines, and in 10 Lacertae there is a line which is 
probably a blend of a strong line of S tv, not heretofore identified, 
with a line of N m1. Nearly all of the sulphur lines measured by 
F. E. Roach’ have been here confirmed. The new identifications of 
sulphur were made with the help of a complete laboratory study of 
this element recently published by M. J. Gilles.? Many new lines 
could be attributed to C 11, the total being 22. Al 1m, which I found 
recently in several B-type stars,’ was confirmed with a total of 15 
lines. Many lines not previously identified could now be attributed 
to Alt. The presence of P 111 was confirmed and several other lines 
were attributed to P 11. Ot is represented by 4, and N 1 by 14 
lines. Ca 1 is represented by the H and K lines. There is a possi- 
bility that two lines which seem to reach maximum intensity in 
6 Cephei are due to Ca 11. Mg 11 is strong throughout the sequence 
of B stars; 7i11 and Fe 1 appear only in the later subdivisions. 
Table II contains one line, \ 4559, which Baxandall has attributed 
to Cr 11, and two lines which he has attributed to Ni 1. 

An attempt has been made to detect other elements which, from 
their excitation and ionization potentials, might be expected to be 
present in B-type stars. There is one line, \ 4186.83, which is per- 
haps identical with \ 4186.24, the strongest line of ionized potassium 
within the range of spectrum observed by me; but the wave-lengths 
do not agree as closely as might be desired. The absence of other 
lines of K 11 (ionization potential 31.6 volts)’ is surprising in view of 
the abundance of this element in the crust of the earth. No evidence 
was found for Bi (ionization potential 24.2 volts) or for Cli 
(ionization potential 22.5 volts). It is also surprising that stellar 

72,191, 1930. 

2 Recherches sur la structure des spectres atomiques du soufre. Paris, 1930. 

3 Astrophysical Journal, 71, 67, 1930. 

4 Ibid., 71, 150, 1930. 

5 T. L. de Bruin, Proceedings of the Royal Academy (Amsterdam), 29, 713, 1920. 
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spectra contain no lines of Ti 11 (ionization potential 27.6 volts),' 
Ti iv (ionization potential 43.1 volts),? or Na 1 (ionization poten- 
tial 47 volts). 

The number of unidentified lines in Table II is not large, but there 
are some which are fairly strong and which appear in many stars. 
Their behavior with respect to type, and consequently with respect 
to temperature, may be ascertained from the intensities in Table II. 
The most conspicuous unidentified line is \ 4419.62. It is observed 
throughout nearly the whole sequence of B stars and reaches maxi- 
mum intensity at or near B2. 


I am greatly indebted to Dr. W. S. Adams for permission to use 
the Coudé spectrograph attached to the 1oo-inch telescope of the 
Mount Wilson Observatory and to Dr. Th. Dunham, Jr., for two 
plates of y Pegasi which have an important part in this investiga- 
tion. Various members of the staff of the Yerkes Observatory have 
co-operated in the observations with the 4o-inch refractor, notably 
Messrs. W. W. Morgan, H. F. Schwede, F. E. Roach, and F. R. 
Sullivan. Acknowledgment is also due, for various suggestions, to 
Professors T. Takamine and Y. Ishida, of Tokyo, and to Dr. R. 
Ruedy, of Toronto. 


YERKES OBSERVATORY 
May 15, 1931 


"H. N. Russell and R. J. Lang, Astrophysical Journal, 66, 13, 1927. 
2 Ibid. 
3 F, H. Newman, Philosophical Magazine, 5, 150, 1928. 


A NUMERICAL METHOD OF DETERMINING THE 
SPACE DENSITY OF STARS' 


By FREDERICK H. SEARES 


ABSTRACT 


The article develops and illustrates a numerical method of solving by approximations 
the fundamental integral equation of stellar statistics for the determination of the 
space-density function. The unknowns are the relative densities d in successive spherical 
shells whose thickness corresponds to a constant increment in the logarithm of the dis- 
tance. The first approximation is obtained by using the fact that the star ratio for mag- 
nitudes m+1 and m, which is known from star counts, bears a simple relation to the 
weighted mean d for the shells which contribute stars of magnitude m. The final approxi- 
mution is made by calculating the effect of changing each d in succession by a standard 
amount, and then combining the series of differential changes, each multiplied by an 
appropriate factor, in such a way as to reproduce the residuals left by the preceding ap- 


proximation. 

No assumption is made as to the form of the density function or the frequency func- 
tions for apparent and absolute magnitudes. The /uminosity function may vary with 
the distance. The illustrative example, solved with a luminosity function which is the 
same at all distances, /eads to a rapid change in density close to the sun that is not sup- 
ported by other evidence. The discordance is discussed in Contribution No. 437. 


I. INTRODUCTION 


The determination of the number of stars per unit volume at 
different distances from the sun depends on the integral equation 


in which w is a constant and N and ¢ are distribution functions for 
apparent (m) and absolute (M) magnitudes, respectively. D, here 
regarded as depending only on the distance p, is the required density 
function. In case space absorption is not negligible, D represents a 
fictitious distribution of stellar density, from which the true distribu- 
tion A may be found, provided that the absorption function is 
known.’ 

Most of the solutions proposed for equation (1) depend on some 
simplifying assumption as to the form of the frequency functions for 


' Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 436. 
2 Mt. Wilson Contr., No. 428; Astrophysical Journal, 74, 91, 1931. 
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N and ¢ or of the density function itself, and none gives a complete 
representation of the observational data. Schwarzschild’s' general 
solution of integral equations of the form of (1) is mathematically 
sufficient, but requires developments in Fourier’s coefficients for V 
and ¢ and leads to extensive calculations if these functions are to be 
accurately represented. Much simpler is his special solution—a 
quadratic exponential in log p for D—which corresponds to the case 
of normal frequency functions for V and ¢. This solution has been 
much used because it permits a simple analytical treatment of many 
problems of stellar distribution. An obvious defect, which implies 
that V and ¢ cannot be exactly represented by Gaussian functions, 
is that close to the sun the calculated densities are much too small, 
and zero at the sun itself. 

In itself, this failure is not very serious, because the density near 
the sun may be found in another way. A greater objection is the 
inflexibility of the method, which requires a Gaussian distribution 
for the star counts. This condition is satisfied, and then but 
approximately, onty when the irregularities in N arising from the 
galactic star clouds and from other stellar aggregations have been 
smoothed out by averaging over many square degrees; but this pro- 
cedure of course also eliminates all irregularity from the calculated 
density function, and the true distribution of density thus escapes us. 

A method which avoids these difficulties begins by resolving equa- 
tion (1), for each value of m, into elements corresponding to moder- 
ate intervals in p (actually, constant intervals in log p are used). 
Since the frequency @ and the volume factor wp’dp are known for 
each element, we thus have, instead of (1), a group of linear simul- 
taneous equations in which the mean densities, D,, D., D;,.... , 
in a series of successive shells, are the unknowns to be determined. 
The practical difficulties of dealing with fifteen or twenty equations 
involving an equal number of unknowns are overcome by a meth- 
od of approximation which leads to a satisfactory representation of 
the star counts. 

It must be noted, however, that the density problem involves an 
inherent uncertainty which apparently cannot be overcome by any 
statistical method of treatment. Because of the enormous range in 


t Astronomische Nachrichten, 185, 81, 1910. 
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the intrinsic brightness of stars, the apparent magnitudes in a 
cluster or star cloud show a range of 15 or 20 magnitudes. The 
sharp narrow maximum in the density function caused by the 
presence of the cloud therefore appears in the star counts as a very 
broad flat maximum. The location of such a cloud is therefore 
always difficult. For a remote cloud the maximum itself may not be 
accessible to observation; the solution then becomes highly sensitive 
to systematic error in the magnitude scale. For a succession of 


TABLE I 
MeEpIAN M CorRESPONDING TO MEDIAN VALUES OF m AND OF LOG p 
MEDIAN m 
MEDIAN Loc p 

0.5 1.5 2.5 3-5 4.5 5-5 

° +1 +2 +3 +4 5 
—2 —I +1 +2 +3 
—3 —2 —I ° +1 +2 
—4 —3 —2 —I +1 

—7 —6 —5 —4 —3 —2 


clouds seen in the same direction the solution is practically indeter- 
minate. For example, the superimposed maxima produced by clouds 
at distances of 1600, 2500, 4000, and 6300 parsecs would be separated 
by intervals of only one magnitude, and without supplementary in- 
formation the resultant curve could not be resolved into its com- 
ponents. The point at which the nearest of a group of clouds begins 
to affect the density function can be located, but what lies beyond 
is likely to remain doubtful. Although the outlook for a complete 
solution of the density problem is not hopeful, the sensitiveness of 
the calculated densities to small changes in the values of V makes it 
the more important that the method used be one which gives the 
closest possible representation of the star counts. 


Il. THE METHOD 


It is advantageous to introduce as unknowns, not the values of D 
itself as suggested above, but the logarithms of the ratios of D for 


™ Cf. Mt. Wilson Contr., No. 347, Introduction; Astrophysical Journal, 67, 125, 1928. 
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successive shells at regularly increasing distances from the sun. For 
these unknowns a first approximation is easily found, and the range 
in their numerical values always remains moderate. For dm in equa- 
tion (1) we use an interval’ of 1 mag., and for dp the corresponding 
interval of o.2 in log p. The median values of M found from 


M=m-— 5 log p (2) 


(Kapteyn’s zero-point for M; p in parsecs) for successive median 
values of m and log p are shown in Table I, which is to be thought 
of as extending downward and to the right as far as necessary. The 
lines of this table correspond to shells having boundaries at log p= 
OD; The volumes of the shells included within the 


angular area w are 
wp'dp=V_s, Vi, (3) 


where the subscripts are 10 times the median log p of the respective 
shells, and where, further, 


V,=oV_, V;=oN,, V,=oV; (log o=0.6) (4) 


If the mean densities of the shells be denoted by D_,, D,, D,; 
. , equation (1) may be replaced by 


5 


or 


N(m) } 


where 


Equations (5) represent the total number of stars in the interval 
m+o.5, the individual terms on the right denoting the contribu- 
tions of the successive shells to this total. For magnitude m+1 we 


have, similarly, 
= +6 Vi © 


« In special cases it may be desirable to use intervals of 0.5 in m and M, and of o.1 in 
log p. The adaptation of the formulae to this case will be obvious. 
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Introduce now the density ratios, d, defined by 
D,=4,.D_,, D,=4d.D,, D;=4d,D, (7) 


where d,, d,,.... refer to the distances log p=o.0, 0.2,..... 
With the aid of (4) and (7), equation (6) then becomes 


N(m+ 1) =a( +) (8) 


and we have the important relation 


oN (m) 9 


This equation states that the star ratio for magnitudes m+1 and m, 
divided by the volume ratio a, is equal to the weighted mean (6,,) of the 
relative densities (d,) for the shells which contribute stars of magnitude 
m. The star counts therefore give the values of log 6,, log 6n4+:, 
fea, coed The problem is to find from the weighted means the values 
of log log d,, . ...... Results for log D_,, log D,, . . . . may then 
be found by a simple summation of the values of log d. 

For the ideal case in which all the distribution functions are 
quadratic exponentials (Schwarzschild’s simplified solution), it can 
be shown that log 6,, equals the value of log d at the distance corre- 
sponding to the mean log p for the stars of magnitude m. Although 
this simple relation does not hold in practice, it suffices to give a 
first approximation for the values of log d,. This we obtain by 
plotting the mean logp for stars of magnitude m, m+1,.... 
against the values of log 6, log bnii,.--.- From the curve thus 
defined, with the arguments log p=o.0, 0.2,. ... , we read the re- 
quired approximations for log d,, log d.,..... 

The mean log p for stars of magnitude m depends, however, on the 
density function, and we therefore need at the outset an approxima- 
tion for D. To illustrate the procedure it is desirable to introduce at 
this point an example. , 

ul. AN EXAMPLE 

For numerical illustration we use the star counts for latitude o° 

given in Table V of Groningen Publication No. 27 and afterward 


4 
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slightly modified by van Rhijn.' Smoothed values of log N(m), 
which refer to the limits m+o.5 and an area of 1 square degree, are 
in the second column of Table II. The results for log 6, found by 
equation (9g) are in the fourth column of the table. Thus, 


log 6,.;= 6.953 —6.411—0.6= —0.058 . 


In addition to these fundamental data, we require values of the 
luminosity function ¢ (/). For these we use the results of Kapteyn 
and van Rhijn,? for which the interval is M+o.5. 

TABLE II 


STAR COUNTS AND VALUES OF LOG 6m 


Os- O—C For Loe 6,, (UNIT =0.001) 
MEDIAN m SERVED | O-C 

Loc N(m)|Loc N(m)| SERVED 

1 Sq. DEG. LOG | (2) (3) (4) | (sa) (5) | Final 
0.5 6.411 +2 |-—o.058| + 3| 3| — 3 —4 —3 —3 
8.027 —2 .069 16} + 3 +2 +3 +3 
4.5 8.558 fe) .074 20|+4/|+1 +2 —I —I 
9.603 .OgI 37 | +14 | + 3 ° ° ° 
.105 42 | +17| +4 ° +1 +1 
0.606 ° .127 43 | +19 | + 4 +1 ° +1 
1.079 ° | — 2 ° +1 +i 
10.5 1.526 +1 .184 32 | +6] — 6 ° —I ° 
1.942 +2 22}/—1] —1I4 —2 —3 —2 
2.079 +2 265 10| — 3] —12 +1 —3 +1 
14.5..........| 3-014 | +3 |—0.271] +32 | +13 5| +2] 


Consider now the stars of magnitude m as expressed by equations 
(5). The contributions of the successive shells have an approxi- 
mately Gaussian distribution. Hence the mean log p for all the stars 
of magnitude m equals the value of log p for maximum frequency. 
For the brighter stars the distances are such that the corresponding 

* Astronomische Nachrichten, 213, 47, 1921. These calculations were begun some years 
azo when no luminosity function for photographic magnitudes was available. It was 
not therefore convenient to use the Mount Wilson star counts. In the meantime van 
Rhijn’s photographic luminosity function, Groningen Publication No. 38, has appeared; 
but, since the present example is only illustrative of the method, the older data are re- 
tained. Applications based on modern data will be made later. 

2 Mt. Wilson Contr., No. 229, Table la; Astrophysical Journal, 55, 245, 1922. 
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. absolute magnitude is about —4.5, and we assume this value to 
represent maximum frequency for all values of m. The values of 
4 TABLE Ila 
i VALUES OF MEAN Loc p 
m (1) (1a) (2) 
“ace 2.4 2.39 2.36 

BOP 4.2 2.72 
TABLE III 
VALUES oF Loc d 
APPROXIMATION* 
Loc p 

: (r) (2) (3) (4) (5) 

—0.058 .056 .O51 .O51 .046 
| .060 .055 .055 .042 
.065 .065 .059 .059 .054 
.069 .070 .062 .062 .072 
q .074 .075 .066 .066 

.082 .083 .O71 .070 .085 
H .OgI .099 .075 .074 .079 
-105 . 130 .084 .077 .059 

.185 . 102 .083 .056 
153 235 .198 . 102 .080 
.218 . 290 .271 .346T 
.245 . 292 . 290 .461 
4 304 .410 364 + .088 
398 . 382 + .256§ 
. 330 .458 .418 . 400 .400 
4 .482 .438 .420 — .420 

—0.358 —0o.506 —0.458 —0.440 —0.440 
: ~* Values of log d for log p<1.0 and> 3.8 are extrapolations. See Figures 1 and 2. 
‘ t Final value is —o0.362. § Final value is +-0.199. 
mean log p corresponding to m=o.5, 1.5,...: are therefore 1.0, 
1.2,.... (Table I) and are to be correlated, term by term, with the 


| 
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data for log 6,, given in the fourth column of Table II. The results, 
which represent the first approximation for the relation between 
log p and log d, appear in the first two columns of Table III and 
as curve 1 in Figure 1. These values of log d give in turn a first 
approximation for D, which, by means of (5), may be used to im- 
prove the values of mean log p. 

For use with (5) it is convenient to form at once the values of DV. 
Final choice of units will be made later; for the present we assume 
both D, and V, (shell having median log p=o0.9) to be unity. Equa- 
tions (4) and (7) then give 


D,V5=1, 

, 
, 


(10) 


The successive values of log DV may therefore be found by adding 
log « =0.6 to each log d in the second column of Table III and form- 
ing the cumulative sums. Since these values of log d are the respec- 
tive values of log 6,,, it happens that this first rough approximation 
for log DV is the original sequence of values for log N(m) modified 
by a constant. Values of log DV for the final approximation, based 
on log d in the last column of Table III, are given in Table VI. 

In using equations (5), Table I is the guide for combining DV 
with ¢. Thus for m=1.5 and log p=1.1, M=~—4. The value of 
D,,V is therefore to be multiplied by ¢ (—4). 

The revision of the values of mean log p does not require a separate 
calculation for each value of m; nor is it necessary to compute each 
term in the right member of (5). Since the relation between m and 
the corresponding mean log p is nearly linear, an interval of three 
or four magnitudes is sufficient. Further, since the distribution of 
is approximately Gaussian, the values of log define a parabola 
whose vertex coincides with maximum frequency. Hence if x,, x., 
x, be three successive values of log ” (preferably near the maximum), 
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and if y, be the value of log p corresponding to x,, the vertex value 
of log p, for the intervals used here, is 


I 
— (%2— 43) (x1) 


Results thus found for six values of m are given in the third 
column of Table IIa. These, together with interpolated data for in- 


log d 
0.0 
—o.1 
—0.2 
—0.3 
la 
>. 
2c 
—0.5 
0.0 1.0 2.0 3-0 4.0 log p 


Fic. 1.—Preliminary approximations for log d 


termediate values of m, are plotted against log6, as curve 1a, 
Figure 1. This curve is the basis for further improvement of log p, 
which leads finally to the last column of Table Ila and to curve 2 of 
Figure 1. The values of logd in the third column of Table III were 
read directly from curve 2 for log p=o.o, 0.2, . . . . and represent 
the results of the second approximation. 

The test of this approximation, as of all others, is afforded by a 
calculation of values of log 6, by means of equations (10), (5), and 
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(9). When the results thus found agree with those derived directly 
from the star counts, the approximation for d, and hence for D, may. 
be adopted as satisfactory. The manner of using (10) has already 
been described. In the case of (5), the calculation is to be extended 
over all shells which yield any appreciable contribution to stars of 
the magnitude in question. As before, Table I is the guide for com- 
bining ¢ and DV. The arrangement of the results is illustrated by 
Table VII, which gives the values of » for the final approxima- 
tion. 

The values of log 6,, found from the second approximation for log d 
are plotted as curve 2c in Figure 1. The corresponding residuals in 
the fifth column of Table II are all positive and show that for the 
third approximation the densities must be increased. The adopted 
corrections for log d are the differences in the ordinates of curves 2 
and 2c; the corrected values in the fourth column of Table III, in 
turn, lead to residuals (3) in Table II. A repetition of the process, 
with the relation between residuals (2) and (3) as a guide, gives the 
results found in Tables II and III under the heading (4). 

The representation is now satisfactory to the sixth magnitude. 
The fourth approximation, however, is no real improvement on the 
third, and some other method of correction must be found which . 
will bring to light the fluctuations obviously present in the values of 
logd. Before this is attempted, a remark may be made concerning the 
revision of calculations involving equations (10) and (5), which may 
often be simplified by noting that the cumulative sums of the correc- 
tions to log d are directly the corrections to log DV. Hence, correc- 
tions to m may be easily computed, and it is frequently advantageous 
to form these corrections instead of completely re-solving equa- 
tions (5). 

IV. APPROXIMATION BY STANDARD VARIATIONS 


The residuals from the fourth approximation are to be reduced by 
further correction of log d to quantities comparable with the accu- 
mulated error of calculation,’ say to two or three thousandths in logé,. 

* This limit, which is far beyond the precision of the data, is chosen to illustrate the 
possibilities of the method. Since residuals (4) of Table II are no great improvement 


over those of the preceding column, the final approximation by the method given in this 
section might advantageously have been started after the third approximation. 
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Further, the final residuals must show an accidental alternation 
of sign. A direct solution is out of the question, first, because of the 
labor involved; and second, because inherent numerical uncertainty 
in the equations would transform small residual errors in the data 
into inadmissible fluctuations in the density function. Various ap- 
proximation methods were tried, none of which would converge. 
Finally a procedure was devised which proved successful and much 
more expeditious than might be supposed. 

In solutions based on the ordinary method of differential correc- 
tion, we form the total differential of the function for simultaneous 
changes in all of the variables. In what follows we use a series of 
partial differentials of the function, each corresponding to a standard 
differential change in a single variable. From these we build up a 
series of changes in the variables, which, made simultaneously, pro- 
duce a series of total differential changes in the function equal to the 
given series of residuals. Thus, for the present problem we make a 
standard change of +0.05 in log d, in succession for each value of 
log p, and calculate in each instance the corresponding changes in 
log 6. The required solution is then found by combining these 
standard residuals, multiplied by appropriately chosen factors, in 
such a way as to reproduce the original residuals. 

Let 


d’=(1+a)d, (12) 


where d and d’ represent original and changed values of the relative 
density. Equations (5) and (10) then give for the corrected number 
of stars of magnitude m 


N'(m) =N(m)+a>-n(m) , 


where the summation begins with the shell whose inner surface 
marks the point at which d has been changed. Similarly, for magni- 
tude m+1 the same change, made at the same point, gives 


. 
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Hence, by (9), the corrected 6,, is 


(13) 
where 


n(m+1); RP = vm) n(m). (14) 


For the standard change 
+A log d,=log(1+a) , (15) 
we have therefore the standard residuals 


+ A log 6,=log (1+ —log(1 +R") , (16) 


where the index m is to be extended over all magnitudes which yield 
appreciable values for the residuals. 

If the standard positive and negative changes in log d are nu- 
merically equal, the respective values of a will not be exactly the 
same. It is simpler, however, to use a constant a and allow the posi- 
tive and negative changes to differ. Thus for the change —o.050, © 
a= —o.109; for a=+o0.109, the corresponding standard positive 
change is +0.045. 

To apply these formulae we begin with the table for m derived 
with the fourth approximation (the equivalent of Table VII) and 
form for each column the cumulative sums from the bottom upward. 
Multiplication of these sums for column m by a/N(m) gives the val- 
ues of R”, where s changes from line to line of the table. For ex- 
ample, values for the line log p= 2.1 refer to the change made at the 
inner boundary of the corresponding shell, that is, at log p= 2.0, for 
which s=20. The remainder of the calculation by equation (16) 
needs no comment. 

Table IV illustrates a portion of the table of residuals (unit = 
0.0001) corresponding to changes of +0.045. Thus if log d for log p 
= 2.4 be increased by this amount and no other change be made, 
the computed values of log 6,, will be increased by the quantities 
given in the first line of Table IV, while residuals (4) of Table II 
will be diminished by the same amounts. 
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The use of these data to build up a solution which will represent 
the star counts depends on the fact that moderate variations in 
the standard changes produce approximately proportional varia- 
tions in the standard residuals. Thus if P times a certain series of 


TABLE IV 


STANDARD RESIDUALS: VALUES OF ALOGédm CORRESPONDING TO A CHANGE OF 
+0.045 IN Loc d at THE DisTaNcE LOG p (UNIT=0.0901) 


m 

1.5 | 2.5 | 3-5 | 4.5 | 5-5 | 6.5 | 7-5 | 8.5 | 9.5 | 10.5 | 11-5 12.5 |13.5|14.5 
4 | 10 21 | 32 | 47 | 61 | 65 | 60 | 49 | 37 | 26 | 17] 8 | 6) 1 
.-.-| § | 10 | 16 | 30] 45 | 55 | 64] 55 | 50 | 41 | 31 | 20] 13 | 6 
2.8. 1] | 14 | 26 | 37 | 48 | 54 | 57] 53 | 47 | 35 | 27 | 18 
3:0. . 1} 2] 3] 11 | 19 | 27 | 37 | 47] 53 | 56 | 53 | 46 | 35 


standard residuals is used, a must be replaced by fa, where f is de- 


fined by 
I + I 
P log ( 08 ) (17) 


The corresponding factor Q for the standard change in log d is then 
given by 
Q log (r+a)=log (1+fa) . (18) 


Expansion of the logarithms gives the approximate relations 


P(+1—4aq) = +f(1 F 4a/q) 
Q(+1-4a) =+f(1F haf) 


where 


These relations show that in general Q differs from P by a quantity 
of the order of a and that for values equal to unity or less the 
difference between P and Q may -be ignored. For cases in which 
large factors are required and for which the proportionality of 


| 

| 
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changes does not hold, the standard residuals may still be used 
by a simple modification of the procedure which will be illustrated 
later. 

The problem is now to find by trial a group of factors such that 
when multiplied into the respective series of standard residuals, 
positive and negative, the sum of the modified series reproduces 
residuals (4) of Table II. The start is made by noting that a positive 
change is required near m= 7.5, and a negative change for m=9.5- 
13.5. Reference to Table IV shows that the maximum effect for 


TABLE V 


FACTORS AND CORRESPONDING CORRECTIONS TO LoG d 
(Fifth Approximation) 


Log p Factor A log d Log p Factor A log d 
ce +o.1 + .004 — 1.5 — 
—0.3 — .O15 0.0 . 000 
— .005 + 3.0 + .123 


these magnitudes is produced by changing log d near log p= 2.6 and 
3.2, respectively. In choosing factors, allowance must be made for 
the fact that overlapping series with opposite signs partially neu- 
tralize each other; further, the curve of corrections must be made 
as smooth as possible by distributing large changes over several 
adjacent values of log d. The criterion of continuity in log d is a 
helpful guide, and a diagram like Figure 2 facilitates the choice of 
the factors. Those used in the present case are given in the second 
column of Table V. The corresponding changes in logd in the 
following column were mostly obtained by multiplying the positive 
and negative factors by the standard changes +0.045 and —0.050, 
respectively. For the cases involving large factors, values of f were 
computed by (17), which were then used with (18) to obtain the 
corrections to log d. The corrected values of log d are under heading 
(5) in Table III. 
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The final residuals to be expected as a result of these changes are 
given by 


Residual (5a) = Residual (4) —2PX Standard Residual. —_(21) 


log d 


+o.r ! —- 


I 


—0.4 


—0.5 
0.0 1.0 2.0 3.0 4.0 log p 


Fic. 2.—Final approximation for log d 


The results, which appear under heading (5a) in Table II, are all 
small and show a proper alternation of sign. 

The real test, however, is the representation of the observed values 
of log 6,, afforded by equations (5) and (9). The residuals from this 
test, (5), Table II, are essentially the same as (57), as they should 
be, except for the last three values. The disagreement here arises 
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from equation (21), which assumes the relation between P and 
Q to be constant, whereas equations (19) and (20) show that it 
varies with m. When the factors are large, the resulting error is not 
inappreciable. An inspection of the calculation suggests the adjust- 
ments at log p= 3.2 and 4.4 given in the footnote to Table III. The 


TABLE VI 
FINAL DENSITIES 
Median Log p Log DV Log D ’ Log Do S—K and oR 
— 2.805 +0.195 —o.781 +0.57 
2.236 0.164 o.812 + .54 
1.671 0.129 0.847 + .50 
1.110 0.090 0.886 + .46 
6.7: —0.553 +0.047 0.929 + .43 
0.000 0.000 0.976 + .37 
+0.554 —o.046 1.022 + .33 
0.088 1.064 + .29 
To Bie 1.657 0.142 1.118 + .23 
2.185 0.214 1.190 + .16 
2.704 0.295 1.271 + .08 
3.219 0.380 1.356 — 
3.740 0.459 1.435 — .08 
4.281 0.518 1.494 — .09 
2 4.825 0.574 1.550 — .04 
5-340 0.654 1.630 + .05 
0.885 1.861 + .04 
5.953 1.248 2.224 — .0 
6.092 1.709 2.685 — .21 
6.279 2.122 3.098 — .26 
6.554 2.447 3-423 — .20 
6.932 2.669 3.645 + .06 
7.620 2.581 3-557 +0.63 
8.410 2.382 3.358 +1.36 
8.619 2.782 3.758 +1.59 
8.799 3.202 4.178 +1.42 
+8.950 — 3.642 —4.618 +2.08 


corresponding values of log DV are in Table VI; m and N(m), in 
Table VII; and the final residuals for log 6,, in the last column of 
Table II. 

An even better test is the representation of the original star 
counts in the second column of Table II. These counts refer to an 
area of 1 sq. deg., but since V represents the volume of a complete 
shell, the values of V(m) in Table VII (unit indicated in last line) 
apply to the whole sky. Further, the units for both D and V have 
arbitrarily been assumed to be the density and the volume of the 
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shell for which median log p=o.9. To make the results comparable, 
the observed counts may also be referred to the whole sky, while the 
values of N(m) in Table VII are multiplied by a factor c which 
transforms them into actual numbers of stars. The comparison itself 
gives log c = 2.522, and finally the residuals in log V(m) shown in the 


TABLE VII 


RELATIVE NUMBERS OF STARS OF MAGNITUDE m AT 
DIFFERENT DISTANCES FROM THE SUN 


MEDIAN m 


MEDIAN 


Loc p 
©.5| 1.5] 2.5] 3-5] 4-5] 6.5] 7.5] 8.5] 


. . [3191 4480] 1509] 4978] 1606) 5009/1 486/41 50] 1081 
I I | 10 | to | | 107] 10 


third column of Table II. The small systematic trend in differences 
reflects the slight tendency toward clustering shown by the positive 
values in the last column of the table. 


V. FINAL DENSITIES—-NUMBER OF STARS PER CUBIC PARSEC 


The cumulative sums of the final values of log d in Table III are 
the adopted values of log D (Table VI), for which, as heretofore, 
D=1 for the shell having median log p=o.9. The density expressed 
in stars per cubic parsec within any shell may be found from the 


| 
20 2 2 
{ +o.1....] §2 7 9 I 
S84) 44) 69] x01 I I 
0.7....| 340] 77] 25] 35) 4) 4 
i 0.9....| 454] 121] 276} 54] 80] 13] 15 2 I 
I.1....] 518] 163] 434] 90] 193] 32] 45} 5] 6 
| 1.3...-}| 509] 187] 588] 157] 357] 70] 116] 16) 20 2 
1.5....] 415] 656] 206] 551] 125] 244] 41| 58 7 I 
1.7....] 277] 140] 603] 221] 186] 423] 82] 137} 19] 23 2 2 
1.9....] 154] 92] 462] 190] 731] 230] 614] 140] 272] 45] 64 8 8 7 
2.1....| 73] 51] 300] 151] 652] 239] 752| 201] 457] 80] 148] 21] 25] 26 2 2 
2.3....| 30] 24] 168] 100] 502] 216] 794] 250] 152] 296] 49} 70} 84 9 8 
2.5....| Io] 10} 84] 58) 346) 175] 276) 867] 232] 527] 103] 171] 242} 29] 30 
3 4| 36] 29] 204] 121| 611] 263] 966} 303] 811] 184] 360] 597] 85] 102 
} 13] 12] 77] 68] 402| 203} 873] 321|1007| 269] 612)1194| 198] 281 
4 3] 28] 23] 158] 04] 474] 204] 750] 236] 279] 463 
5 5} 27| 163} 82] 353] 130] 407|r089| 248] 483 
I 7 5| 38} 22] 113} 49] 170] 561] 150} 341 
| N(m).. 
Unit... 
q 
| 
j 
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corresponding line of Table VII. For any such line the horizontal 
sum is 


S= n(m)=DVs , s= 6(M) (22) 


where m, and m, are the apparent magnitudes within which n 
has appreciable values and M, and M, the corresponding absolute 
magnitudes. To obtain the total number of stars in the shell, S 
must be multiplied by the constant c of the preceding section; fur- 
ther, V is to be expressed in cubic parsecs. The required density 
then becomes 

S 


DAo= 


(23) 
where A, is the number of stars per cubic parsec in the shell for 
which D is unity. The constant A, may be computed from the data 
for any shell, or, as a check, from two or more shells. Thus for log p= 
0.9, 1.9, and 2.9, log A,=9.026, 9.024, and 9.025, respectively. 

The corresponding values of log D, in the fourth column of Table 
VI may be compared with the results of Kapteyn and van Rhijn 
(log A). From about too parsecs outward their densities were de- . 
rived from exactly the same data as used here, on the assumption 
that A has the form of a quadratic exponential in log p. For the inter- 
val 100-2000 parsecs, the values of log D,/A in the fifth column of 
Table VI correspond to an average difference of about 12 per cent. 
At greater distances the difference in the two solutions is much 
larger. The uncertainty in the data is such, however, that probably 
there is little choice between the two series of results; and the point 
made is merely that the present solution closely represents the 
initial data. Near the sun, log A depends on different data, treated 
by a different method. The practical constancy of A in this region 
is confirmed by a later investigation by van Rhijn? and cannot be 
seriously questioned. The significance of the discordant values of 

« Mt. Wilson Contr., No. 229, Table Ib, second column; Astrophysical Journal, 55, 
246, 1922. 

2 Groningen Publication, No. 38, p. 15, 1925. 
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log D, in the vicinity of the sun is discussed in the following Contribu- 
tion; here it need only be noted that D, represents the star counts. 


VI. FINAL REMARKS 


In conclusion, it may be noted again that the method imposes no 
restriction on the form of the distribution functions for apparent and 
absolute magnitudes. Further, allowance for any dependence of 
the luminosity function on distance is easily made by using in equa- 
tions (5) the frequencies for M which are appropriate to the re- 
spective shells. The illustrative example has been treated with a 
constant ¢ because hitherto we have had no reason for doing other- 
wise. The rapid change in the resulting densities in the neighbor- 
hood of the sun indicates, however, that this function is not every- 
where the same. See Contribution No. 437. 

The necessity for extrapolating the curves for log d suggests a 
certain arbitrariness in the solution. The results do depend on the 
extrapolation, but the analytical methods likewise involve an extra- 
polation in that the function which represents the star counts is 
tacitly assumed to represent fainter stars not included in the counts. 
Table VII illustrates the difficulty. Magnitudes 12-16, for ex- 
ample, include a certain percentage of stars which are more than 
10,000 parsecs distant. The distribution of these stars depends on 
the density at great distances, which for its precise determination 
would require counts of very faint stars. Since these counts are not 
available, extrapolations of counts of the brighter stars must be 
used. Table VII indicates, however, that up to 10,000 parsecs the 
uncertainty arising from the extrapolation is not very serious. The 
total number of stars in any column is determined by observation; 
the percentage beyond 10,000 parsecs is not large. Hence some un- 
certainty in this percentage does not greatly affect the distribution 
of the nearer stars of the same apparent magnitude. 

Although the initial star counts are well represented, the calcu- 
lated densities are not strictly correct because no allowance has 
been made for certain small statistical errors which enter at various 
points in the calculation. For example, the values of M in Table I 
correspond to the median values of log p and of m within certain 
intervals and do not exactly equal the mean values of M within these 
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same intervals. Again, the mean density of a shell does not occur at 
a distance corresponding to median log p. Correction is easily made 
whenever it is worth while to do so. 

Finally, it may be noted that statistical relations between ap- 
parent and absolute magnitudes, distance, and parallax—for ex- 
ample, the mean parallaxes of stars of magnitude m—may easily 
be formed from the array of numbers in Table VII. 

I am greatly indebted to Miss Joyner of the Observatory staff for 
much assistance with the computations required for the foregoing 
discussion. 

CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
August 1931 
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REVIEWS 


Wavelength Tables for Spectrum Analysis. Second edition. By F. 
TwyMan and D. M. Situ. London: Adam Hilger, Ltd., 1931. 
Pp. xii+180. 14s. od., post free. 

Those who have used the first edition of the Wavelength Tables will find 
in this edition a much improved book. Besides a revision of the wave- 
lengths which were used merely for the identification of an element, the 
authors have added information regarding, and methods of making, 
quantitative analyses with the aid of spectroscopy. As an introduction 
for the practical worker to the theory underlying modern methods a sec- 
tion on “The Various Types of Spectrum” by Professor E. N. da C. 
Andrade has been inserted. The new edition considers analysis by the 
arc and flame methods and has a table of the ‘Sensitive Arc Lines of 
50 Elements” which lists the lines according to wave-length for each ele- 
ment, gives the intensities and the relative sensitivities of the lines. An- 
other section deals with the characteristics and use in analysis of the 
sensitive lines in the arc and spark spectra. This section contains a classi- 
fication of the ultimate lines of 61 elements, the table containing the wave- 
length, the intensities in arc and spark, the series relations, the quantum 
number Ak, the authorities, the excitation potentials, and the tempera- 
ture classification of A. S. King. 

Part IV of the first edition, ‘Wavelengths for Use in Determination of 
Stellar Radial Velocities,” has been omitted. This is quite justified, since 
the book was prepared for the use of chemists and metallurgists in spec- 
trum analysis. The astronomers have access to the tables in other publica- 
tions; and besides, it is an opportune time for a book to be written on 
stellar spectra, collecting together a revision of tables of wave-lengths 
suitable for the determination of radial velocities, the identifications of 
lines, and a discussion of the variations of the lines with the physical con- 


ditions in the stellar atmosphere. 
C. T. ELvey 
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